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ABSTRACT 



Context. During the last decade, the FORS1 instrument of the ESO Very Large Telescope has been extensively used to study stellar 
magnetism. A number of interesting discoveries of magnetic fields in several classes of stars have been announced, many of which 
obtained at a ~ 3 a level; some of the discoveries are confirmed by measurements obtained with other instruments, some are not. 
Aims. We investigate the reasons for the discrepancies between the results obtained with FORS1 and those obtained with other in- 
struments. 

Methods. Using the ESO FORS pipeline, we have developed a semi-automatic procedure for magnetic field determination. We have 
applied this procedure to the full content of circular spectropolarimetric measurements of the FORS 1 archive (except for most of the 
observations obtained in multi-object spectropolarimetric mode). We have devised and applied a number of consistency checks to our 
field determinations, and we have compared our results to those previously published in the literature. 

Results. We find that for high signal-to-noise ratio measurements, photon noise does not account for the full error bars. We discuss 
how field measurements depend on the specific algorithm adopted for data reduction, and we show that very small instrument flexures, 
negligible in most of the instrument applications, may be responsible for some spurious field detections in the null profiles. Finally, 
we find that we are unable to reproduce some results previously published in the literature. Consequently, we do not confirm some 
important discoveries of magnetic fields obtained with FORS1 and reported in previous publications. 

Conclusions. Our revised field measurements show that there is no contradiction between the results obtained with the low-resolution 
spectropolarimeter FORS1 and those obtained with high-resolution spectropolarimeters. FORS1 is an instrument capable of perform- 
ing reliable magnetic field measurements, provided that the various source of uncertainties are properly taken into account. 

Key words. Magnetic field - Polarization - Astronomical instrumentation, methods and techniques - Stars: magnetic field 



1. Introduction 

During a full decade of operations, the FORS1 instrument of 
the ESO Very Large Telescope has collected a large number of 
magnetic field measurements in various kinds of stars. Together 
with the ESPaDOnS instrument of the Canada-France-Hawaii 
Telescope, and with the MuSiCoS and NARVAL instruments of 
the 2 m Telescope Bernard Lyot of the Pic-du-Midi Observatory, 
FORS 1 has been one of the workhorse instruments for the ob- 
servational studies of stellar magnetism. 

Several important detections obtained with FORS1 have 
led to the conclusion that magnetic fields are quite com- 
mon in a variety of stars across the Hertzsprung-Russell di- 
agra m, including for in stance central star s of planetary neb - 
ulae dJordanet all 12005), hot subdwar fs dO' Toole et all 120051) . 
jSCephei and slowly puls ating B stars i HubrigetalJ 2009a ). B 
stars with emission lines (jHubrig et all 2009bl) . and normal O- 
type stars dHubrig et alj|2008bl) . 

In practice, a close inspection to the published results shows 
a number of problems: 

i) Inconsistencies between field measurements obtained with 
FORS 1 and fi eld measurements ob tained with other instruments. 
For instance, iHubrig et alJ d2004al) reported the discovery of a 



magnetic field in the Herbig Ae/Be star HD 139614, while re- 
peated ESPaDOnS measurements fail ed to confir m the magnetic 
nature of that star (Wad e et al] 120051) . Similarly. Isilvester et alJ 
(2009) failed to confirm the presence of a magnetic field in sev- 
eral of the Cep and SPB stars observed by Hubrig et al. (2006a) 
and lHubrig etail d2009al) . 

ii) Inconsistencies between the analysis of the sa me FORS 
datase t performed by different authors. For instance, iMcSwainl 
(2008) observed the normal B stars NGC 3766 MG 111 and 
NGC 3766 MG 176, and the Be star NGC 3766 MG 200, and 
reported no field detection. Using the same data, Hubri get alJ 
(2009b) reported new field detections for all three stars. 

Hi) Inconsistencies between the analysis of the same FORS 
dataset perfor med by the same a uthors but at different epochs. 
For instance, IWade et all d2005l) reported a possible detection 
in the young Ap star HD 72 106 A, which was not confirmed 
by the later analysis of th e same data performed by the same 
group dWade et al.l [2007ah . Note that the magnetic nature of 
that star was est ablished with data independently obtained with 
ESPaDOnS (see lFolsom etal]l2008l) . 

iv) Inconsistencies between field measurements obtained from 
different subsets of an observing series of frames. Magnetic 
field measurements are often obtained by combining a number 



1 



S. Bagnulo et al.: Magnetic field measurements and their uncertainties: the FORS1 legacy 



of pairs of frames obtained at two different position angles of 
the retarder waveplate. This redundancy is mainly motivated by 
the need to reach a very high signal-to-noise ratio. In some rare 
cases, a magnetic field may be firmly detected in a pair of frames, 
but not in the remaining pairs. This is for i nstance the case of th e 
measurements of HD 139614 reported bv lHubrig etal . (2004a), 
where the magnetic field is detected only in a subset of frames, 
and in a couple of H Balmer lines. 

v) Finally, there are some global inconsistencies of the full FORS 
dataset, revealed for instance by the high incidence of field de- 
tection in the null profiles; this kind of problem was not previ- 
ously reported in the literature, and will be discussed in Sect. 14. II 

Issues ii) and Hi) point to possible glitches in the data reduc- 
tion method, while issues i), iv), and v) might be the symptom of 
a wider, possibly instrumental, problem. 

T he release of the FOR S pipeline for spectropolarimetric 
data (Izz o et al. I l20loll20ri gave us the opportunity to develop 
an accurate and efficient reduction method. The FORS pipeline 
is a software tool specifically designed for that instrument, which 
makes it easier to handle some characteristics that are specific to 
the data obtained with FORS 1 . Furthermore, the FORS pipeline 
allows a high degree of automatization in the data reduction pro- 
cess. Complemented with a suite of software tools for data pre- 
processing (for frame classification and quality check) and data 
post-processing (for the computation of the magnetic field), we 
were able to build up a nearly automatic tool for data analysis 
which allowed us to treat the entire FORS archive data in a ho- 
mogeneous way. Using our tool suite, the difference in terms of 
effort required to perform the analysis of a single series of obser- 
vations compared to performing the analysis of the entire archive 
of FORS 1 data consists mainly in the amount of the necessary 
disk space and CPU time. 

Compared to the reduction of raw data coming from individ- 
ual observing runs, the mass-production of reduced spectra of- 
fers the possibility to perform a quality check of the final prod- 
ucts on a very large scale. The aim of this paper is to present: 
(1) a preliminary discussion of the methods, (2) the results of 
our quality checks, and (3) to caution the reader about the lack 
of robustness of some of the results previously published in the 
literature. A deeper analysis and a comprehensive and homoge- 
neous catalogue of FORS 1 magnetic field measurements will be 
published in a forthcoming paper. 

2. Instrument and instrument settings 

FORSl ( Appe nzeller et al.ll 1998b is a multi-purpose instrument 
capable of doing imaging and low-resolution spectroscopy in the 
optical, equipped with polarimetric optics. It was attached at the 
Cassegrain focus of one of the 8 m units of the ESO Very Large 
Telescope (VLT) of the Paranal Observatory from the begin- 
ning of operations in 1999 until instrument decommissioning in 
March 2009. The FORSl polarimetric optics were subsequently 
moved to the twin instrument FORS2. 

Two detectors have been used in the FORSl instrument: a 
2k x 2k Tektronix CCD (TK2048), (pixel size 24/mi x24//m), 
and a mosaic composed of two 2k x 4k E2V CCDs (pixel size 
of 15 /imxl5 yimi). TK2048 was used from the beginning of op- 
erations to March 2007 (P65 to P78), E2V from April 2007 to 
the decommissioning of the instrument in March 2009 (P79 to 
P82). F or more details about the two detectors, see lSzeifert et alJ 
d2007h . 

Most of the observing programmes used the so-called "fast- 
mode". In this mode the target is just one star in the centre of 



Table 1. Summary of the characteristic of the grisms+CCD 
setting most commonly employed for magnetic field measure- 
ments. 



Grism 


CCD 


Wavelength 


Dispersion 


Resolution 






range (A) 


ApixeL 1 


(l"slit width) 


600 B 


TK2048 


3480-5900 


1.2 


830 


600 B 


EEV 


3400-6100 


0.75 


810 


1200 g 1 


TK2048 


4290-5470 


0.58 


1725 


1200B 


EEV 


3660-5110 


0.36 


1540 


600 R 2 


TK2048 


5250-7415 


1.2 


1230 



1. We note that grism 1200 g was often used setting the slit close to the 
right edge of the instrument field of view. For that special setting, the 
observed wavelength interval was 3840 - 4970 A. 

2. Grism 600 R was used together with order separation filter GG435. 



the instrument field of view. A small number of observing pro- 
grammes used the multi-object mode, which allows the simulta- 
neous recording of the polarized spectra for up t o nine different 
stars within a ~ 7' X 7' field of view (see e.g. Bagnulo et al. 
2006). In multi-object mode, the FORS pipeline is not yet robust 
enough to correctly associate the parallel and the perpendicular 
beams of the same target, especially when the observations are 
performed with the grism 600 B. The reason is due to the pres- 
ence of a number of ghost images (probably coming from the 
Longitudinal Atmospheric Dispersion Corrector), which are in- 
terpreted by the pipeline as stellar spectra. The solution to this 
problem would be to force the pipeline to extract only user se- 
lected spectra, an option that has not been implemented yet. 
Therefore we did not re-reduce the large majority of the obser- 
vations performed in multi-object mode. 

Most of the observations were performed with grism 600 B, 
some with grism 1200 B, and a small number with grism 1200 g 
and 600 R (see Table [TJ. In most of the cases, the slit width was 
set to 0.4" or 0.5", for a typical spectral resolution of about 1600 
- 2300 (with grism 600 B or 600 R) or 3000 - 3400 (with grism 
1200 g or 1200B). 



3. Data reduction 

In the following, we will adopt the same formalism used in 
Bag nulo et alJ d2009b : /" and f x are the fluxes in the parallel 
and in the perpendicular beam of the polarisation analyser, re- 
spectively, Py — V/I is the circular polarization normalised to 
the intensit y, and Ny is the nul l profile, a quantity that was in- 
troduced bv lDonati et al.l dl997l) . and that is representative of the 
noise of Py, 

We always obtained the Py profiles from a series of one or 
more pairs of exposures. Each pair of exposures is composed 
of two frames obtained with the retarder waveplate at position 
angles separated by 90°. For most of the observing series, it 
was also possible to calculate the null profile. For those cases 
in which the number of pairs of exposures TV was odd and > 3, 
the null profile was obtained omitting the last pair of exposures. 
Obviously, with just one pair of exposure, the null profile could 
not be calculated. 

The extracted fluxes were combined to obtain the Py and Ny 
profiles using the formulas of the difference methods given in 
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Eqs. (A2) and (A7) of iBagnulo etail {2009), which for conve- 
nience we reproduce below: 



N 

2N 

7=1 
N 



Nv = 4 ZC-D " 1 ' 



(1) 



where AA = 90° and a } belongs to the set {-45°, 135°}Q 
However, in this work we adopt an important modification as 
described below. 

3. 1 . Rectifying Stokes profiles 

In some cases we found the Py profile clearly offset from zero, 
even when no circular polarization of the continuum is expected. 
This happened for instance in Herbig Ae/Be stars observed by 
Wad e et al.l (l2007ah . and in several other cases in the course of 
the present work. 

A possible explanation is cross-talk from linea r to circular 
polarization, as discussed by Ba gnulo et al.l d2009l) . Cross-talk 
may be a problem only if the observed source is linearly po- 
larised, and is especially significant for spectra acquired with a 
slitlet close to the edge of the instrument field of view (as in 
some observations obtained in multi-object mode). However, we 
often found slight but noticeable offsets also in FORS data for 
stars that are not linearly polarised, and that were observed in 
the centre of the field of view. These offsets may be explained if 
the ratio between the transmission functions in the perpendicu- 
lar beam h x , and the transmission function in the parallel beam, 
/z", is not constant as the retarder waveplate is set to the different 
position angles. 

Inspection of the null profile may help to discriminate be- 
tween the two cases. If Py is offset from zero due to cross-talk 
from linear polarization (or in fact because the source is intrinsi- 
cally circularly polarized in the continuum), the null profile will 
still be oscillating about zero, while an offset introduced by a non 
constant ratio /z x /7z" will still appear (somewhat scaled down) in 
the null profile. This can be seen analytically as follows. We first 
recall that within the framework of the difference method, the 
ratio h = h x /h" is developed to first order: 



h 



1 +6h 



(2) 



We then assume tha t the term 6h that appears in Eqs. (34) of 
Bagn ulo et al.l (120091) depends on the position angle cyq of the 
retarder waveplate, and we repeat the computations that lead to 
their Eqs. (37). We obtain 



afreet) 



(red) 



— y 

2N H 
J=i 
N 



n: 



where 



a 



( f-r ) (f-r) 

\M% I ' •/ },, 



+AA 



(3) 



2[dh(aj+AA)-Sh(Oj)] 



J 4+2[Sh(aj)+Sh(a j +AA)]+Sh(a J )Sh(a j +AA) 
j ~ 4+2[5h(a j )+5h{a j +hA)\+6h{a 1 )Sh{a j +hA) 



(4) 



1 We note that instead of setting the A/ 4 retarder waveplate to all four 
possible angles, most of the observers preferred to use only the angles 
-45° and +45°. 



and where we have introduced the new symbols P^ ect) and A^ ect) 
to denote the Py and Ny profiles rectified. 

Equation 0) shows that if 5h(aj)5h(aj + AA) <K 1, then 

bj =* 1, and Pf ct) and A^ rect) can be simply obtained by sub- 
tracting from the Py and Ny profiles obtained from Eqs. (0 two 
smooth curves that interpolate the continuum of those profiles. 
These smooth curves may be obtained fitting a polynomial to 
Py and to Ny, or using a Fourier filter. Following this simple 
method, Py and Ny are rectified independently, and Ny loses 
part of its diagnostic content, because it is forced to oscillate 
about zero in a way that is totally independent of Py. An alter- 
native method consists of estimating the transmission functions 
Sh by smoothing the ratios between fluxes on the parallel and 
perpendicular beams, then obtaining the explicit expressions for 
the individual aj and bj coefficient, and finally using Eqs. 01 



to compute the p^ 1 ' and N v profiles in a consistent way. In 
practice, we found that the two alternative methods lead always 
to very similar results. 



3.2. Clipping with the null profile 

Both Py and Ny profiles show occasional spikes that occur in 
the same wavelength bin. Most of these spikes are probably pro- 
duced by cosmic ray hits, and, if not removed, may lead to spu- 
rious detections of magnetic fields, or at least decrease the pre- 
cision of its determination. A possible remedy is to exploit the 
null profiles to clip the Py profiles by discarding from the com- 
putation of the (B z ) values those points for which the (rectified) 
Ny value depart from zero by more than k <x, where k is a con- 
stant typically ~ 3. The clipping may be applied to the individual 
deviating points, or may also remove the adjacent points. 



3.3. Magnetic fields determinations 

FORS longitudinal magnetic field measurements are obtained 
using the relationship 

, dl(A) 

1 (B z ) , 



V(A) 



Jeff <-Z 



C z A' 



AA 



(5) 



where V(A) and 1(A) are the Stokes V and / profiles of a spectral 
line, g e ff is the effective Lande factor, and 

e 



C z = 



4mn e c 2 



(^4.67x KT^A- 1 G" 1 ) 



(6) 



where e is the electron charge, m e the electron mass, and c the 
spe ed of light. 

IBagnulo et al.l J2002) proposed to use a least-squares tech- 
nique to derive the longitudinal field via Eq. ©, by minimising 
the expression 



X 



2 v (yi-(B z ) Xi -b) 2 



-Z 



(7) 



where, for each spectral point i, y,- = Py(Ai), Xj = 
-geff CzAj(l/F xdF/dA)j, F, is the flux measured in the spectral 
bin at Aj, and b is a constant introduced to account for possible 
spurious residual polarization in the continuum. Note that the 
spurious polarization is assumed constant in wavelength, which 
may not be true if we adopt Py from Eqs. (Q]), but it is a real- 
istic assumption if we use P ( ^ ect) of Eq. (0, for which we must 
retrieve b — within its error bar. 

Equation (0 is subject to several important limitations, (i) It 
is valid only in the limit of fields weak enough that Zeeman split- 
ting is small compared to the local spectral line width, typically 
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■a 



x 



.M/mUiaJU.I.Iaa.. 
pvv|l|-wyy\/|"'V"V Tvy 



iaaA kA AA«> ^A^AaA* AA A Aa^ I/Ua Aaa_. A iA a. a_ a..j.a. a. a . jAa A a 



_l 1_ 



_l I 



c 

3 

.£> 

■5- 

a. 



0.5 



A -0.5 



0.5 



-0.5 




3800 



4000 




4200 4400 
WAVELENGTH (A) 



0.5 



4600 



4800 



-0.5 



-lxlO" 8 10" a 
-4.67 10" 13 X 2 (1/F) (dF/dX) (G" 1 ) 




-lxlO" 6 lO- 6 
-4.67 10" 13 X 2 (1/F) (dF/dX) (G->) 



Fig. 1. The polarised spectrum of HD 94660 observed on 2007-03-22 with grism 600B. The top panel shows the normalised 
derivative of the observed flux, and the panel immediately below shows the observed flux F (black solid line, in arbitrary units, 
and not corrected for the instrument response), the Py profile (red solid line centred about 0), and the null profile (blue solid line, 
offset by -0.75 % for display purpose). The null profile is expected to be centred about zero and scattered according to a gaussian 
with cr given by the Py error bars. The Py error bars are represented with light blue bars centred about -0.75 %. The slope of the 
interpolating lines in the bottom panels gives the mean longitudinal field from Py (left bottom panel) and from the null profile 
(right bottom panel) both calculated using the H Balmer lines only. The corresponding (B z ) and (N z ) values are -1885 + 90 G and 
-192 + 62 G, respectively. 



with a 5 1 kG strength in the case of optical metal spectral lines, 
or 55 10 kG for Balmer lines, (ii) It applies to isolated, unblended 
lines, (iii) The value of g e ff varies by a significant amount from 
line to line; the use of an average value means that for many 
lines, the actual circular polarization will vary from the average 
one by u p to 25%. The im pact of these limitations was dis- 
cussed by Bag nulo et alJ (120021) : they are not important if one is 
simply interested in field detection, or to asses the magnitude of 



the longitudinanal field, but they should be kept in mind when- 
ever FORS measurements are used for modelling purposes. 

When we derive (B z ) by minimizing the x 1 of Eq. (0, we 
make the implicit assumption 



1 dl(A) 
7(1) cU 



1 dF 
F dl 



(8) 
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(9) 



which is justified provided that 

| J_^,» | J_^| 
'/(/I) dA T(/l) cU 

where T(A) is a function that accounts for the stellar contin- 
uum, the transmission functions of the interstellar medium, of 
the Earth atmosphere, of the telescope, and of the instrument. 
The validity of this approximation can be empirically evalu- 
ated case by case by verifying that the impact of the contin- 
uum rectification to the normalised derivative of Stokes / is 
negligible. Practically, a quick check can be performed even 
without continuum normalization. Since the function T(A) is 
much smoother than 1(A), the contributions due to the terms 
[l/I(A) x dI(A)/dA] and [1/T(A) x dT(A)/dA] to the quan- 
tity [1/F x dF/dA]i can be visually disentangled simply by 
plotting the function [I IF x dF/d/l], versus A over a zero line. 
The top panel of Fig.Q]shows that the broad contribution due to 
T(A) is totally negligible compared to the fine structure due to 
1(A), and that [1/F x dF/d/l], practically oscillates about zero. 
For hot stars observed with grism 1200 B, T(A) is very steep, and 
[I IT (A) x dT(A)/dA] may depart from zero more significantly 
than what displayed in Fig. [TJ yet the impact of the continuum 
normalisation on the (B z ) determination is still within the error 
bars due to photon-noise. 

Finally, we note that for quality check purpose, it is useful 
to measure the mean longitudinal field by using the null profiles 
instead of the Py profiles. This null field value (N z ) should be 
found consistent with zero within its error bars. 

Figure Q] illustrates the technique discussed in this Section, 
showing the results obtained for the FORS observations of the 
well known magnetic star HD 94660. 

3.4. Error bars 



Bagnu lo et al ] (l2002h have obtained the error bar of the longitu- 
dinal magnetic field from the inverse of the^- 2 matrix. This way, 
the "internal" error bar S(b z ) was obtained by propagating the er- 
rors cr, of the Py profiles onto the analytical expression of (B z ) 
that is obtained by minimising the x 2 expression of Eq. (0: 



\B Z ) 



(10) 



where is a sum extended to all i spectral bins used in Eq. (0. 
An analogous expression holds for the error on the null field 
(N z ), S(n_), which turns to be numerically identical to S(b z ), since 
erf are identical for the Py and Ny profiles. 

Since there are numerous indications that the error bars of 
FORS 1 field measurements are underestimated, in this work we 
introduce an "external" error bar <T(b.) that takes into account the 
actual scattering about the interpolating line, i.e.: 



- S (B Z ) yf; 



X ■ lv 

™ min' 



(ID 



where x\m is the minimum value of the x 2 of Eq. (0 and v is 
the number of degrees of freedom of the system, i.e., the num- 
ber of spectral points minus two. Strictly speaking, the use of 
Eq. (fTTT i is not theoretically justified. It is equivalent to the hy- 
pothesis that Eq. (0 is correct, but that the errors cr, on Py are 
all underestimated by a constant value (y 2 . /v) 1 ^ 2 . 

J va man ' ' 

The reality is that, as discussed above, Eq. ([5]) is not rigor- 
ously valid, especially for strong fields. Departures from the be- 
haviour predicted by Eq. (0 mean that points in the regression 



(as shown in the lower panels of Fig. will deviate from the 
interpolating line by more than the amount predicted by pho- 
ton statistics, and the reduced x 2 values will be > 1 . For stars 
with weak or no magnetic field, or if the field is computed from 
the null profiles, the scatter in the positions of individual points 
around the regression line is dominated by noise rather than by 
the approximations implicit in Eq. (0. In this case the reduced 
X 2 is expected to be very close to 1, and the "internal" and "ex- 
ternal" errors tend to be quite similar. The advantage of system- 
atically using Eq. (fTTI) is to automatically take into account a few 
cases where even if the field is weak or non-existent, the reduced 
X 2 is > 1 . 

To check the validity of our error estimates we have per- 
formed some Monte Carlo simulations. All pixel values of the 
raw frames of an observing series were scattered, using a ran- 
dom number generator, according to a Gaussian distribution with 
cr given by the square root of their electron counts. This altered 
dataset was reduced, and an estimate of the (B z ) and (N z ) val- 
ues was obtained. The same procedure was repeated 5000 times 
to obtain a distribution for both (B z ) and (N z ) values. Then, the 
centres of these distributions (B Z ) MC and (N Z ) MC , and their stan- 
dard deviations crjl c and cr^ c were compared to the (N~), 
<T(B-), V{n z ), S(b z ) (= S(n >) values obtained from the observing se- 
ries not altered by the Monte Carlo simulation. Results for three 
observing series are shown in Fig. [2] For the (large) (B z ) values 
of HD 94660 we found ct^b.) > S(b z ) and <T(b„) > f° r a H the 
remaining cases the error bars coming from the three methods 
are similar. 



3.5. Additional details about data reduction 

In this Section we briefly discuss some issues related to data 
reduction which influence the estimates of the mean longitudinal 
magnetic field and their uncertainties, i.e., i) flat-field correction; 
ii) the algorithm for spectrum extraction; ii) spectral rebinning; 
iv) rectification of the Py and Ny profiles; and v) the choice of 
the Lande factor and of the spectral region used to determine the 
magnetic field. 

This (non-exhaustive) list gives an idea of the multiple 
choices that are either implicitly or expliciting taken during the 
process of data reduction. The following discussion helps to un- 
derstand how much the final results may depend on the actual 
algorithm adopted for data reduction. 

i) While flat-fielding is not n eeded to measure Stokes profiles 
(see, e.g jBagnulo et all2 009), pixel-to-pixel variation may have 
an imp act on the evaluation o f the x, term of Eq. (0, as pointed 
out by iRivinius et alJ ((2010). However, in most of the cases, 
the signal-to-noise ratio of the combined flat field frames is not 
much higher than that of the science frames. Furthermore, cali- 
bration frames show internal reflections (especially in the bluest 
regions), but with different patterns than seen in science frames. 
In some cases, it is clear that internal reflections observed in the 
flat-field affect the (B z ) determination, e.g., the difference be- 
tween cases (a) and (b) of Table|2]described later in this Section 
is entirely due to spurious features of the master screen flat. In 
summary, it is not possible to conclude that flat-fielding cor- 
rection improves the longitudinal field determination, at least 
with the CCDs that were used in FORS1. We found however 
that the difference between magnetic field values obtained with 
and without flat-fielding corrections were generally well within 
photon-noise error bars. 

ii) The FORS pipeline dlzzo et al.l l201 1) offers tw o different op - 
tions for spectral extraction: optimal extraction dHornd [l986). 
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Fig. 2. Results obtained from the analysis of the H Balmer lines (top panels), of the metal lines (mid panels), and of the full 
spectrum (bottom panels) for three stars, using standard error theory and Montecarlo simulations. The (red) thick solid lines are 
Gaussian curves centred about the various field values and with <x given by Eq. (TTTb . The blue lines show the distributions obtained 
with Monte Carlo simulations as explained in the text. Note that HD 94660 has a large field, HD 96441 has no detected field, and 
HD 171 184 has a weak but probably real field. 



and aperture extraction. Using the former method, each pixel 
in the spectra extraction is weighted according to the fraction 
of the flux which is expected assuming a model of the spatial 
profile across the dispersion. The method substantially increases 
the signal-to-noise ratio when the noise is dominated by back- 
ground, and automatically accounts for cosmic-rays, but may be 
less suitable than "average extraction" when the noise is deter- 
mined mainly by Poisson statistics. We found that in general, 
the polarized spectra obtained with optimal extraction have nu- 
merous spurious spikes; an example is shown in Figure [3] with 
the detail of the polarized spectrum of HD 94660 reduced with 
average extraction (blue thick lines) and with optimal extraction 
(red thin lines). This suggests that the aperture extraction should 
be preferred to optimal extraction. Note that FORS pipeline as- 
sumes a certain default aperture which may not be the optimal 
one for all cases. For our data re-reduction we generally adopted 
a 12 pixel aperture, using a larger one for a few cases of very bad 
seeing. 

Hi) Generally, the wavelength bin size of the calibrated spec- 
tra is set equal to the dispersion (measured as A per pixel) per- 
taining to each given instrument setting, but there is no obvi- 
ous reason why one should not experiment with an interpolation 



over a more refined wavelength grid, or with a rebinning to a 
larger bin size to increase the signal-to-noise ratio of each point. 
Surprisingly, the final (B z ) value depends on the adopted choice 
of the bin size, in a way that deserves further investigation. 

iv) The choice whether to use Balmer lines only or the full 
spectrum, and the choice for the Lande factor that best represents 
the average Zeeman pattern of the spectral lines may be different 
from work to work, and therefore lead to different numerical val- 
ues for (B z ). However, since these choices affect both the field 
value and its error bar in the same way, they will not change a 
detection into a non-detection nor vice-versa. For more details 
on how field results change when Eq. (0 is applied to Balmer 
lines onl y, or only to He+me tal fi nes, or to the full sp ectrum, 
see, e.g., [Bagnulo et al. (2002) and Bagnulo et al ] (l2006h . 

v) In some cases, using the rectified profiles of Eq. (f3j) instead 
of the non-rectified Py profiles of Eq. (Q]) leads to significant 
changes in the field values, and may even turn a field detection 
into a non-detection, but in most of the cases, field changes are 
within the error bars. 

Table [2] shows the field values calculated for three stars us- 
ing various algorithms. It is interesting to note that a field mea- 
surement of the Ap star HD 171 184 could be found to be either 
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Fig. 3. P v and N v profiles of HD 94660 obtained with aver- 
age extraction (red thick lines) an optimal extraction (blue thin 
lines). The null profiles are offset by -1 % for display purpose, 
and overplotted to the statistical error bars of the Py and Ny pro- 
files. 



as a non detection, a marginal detection, or even a Act detec- 
tion, according to the selected reduction algorithm (since the 
star is chemically peculiar, probably the star has a weak mag- 
netic field). As a further example of the impact of a choice made 
during data reduction, Fig. [4] shows the results the distribution 
of the difference between (B z ) values obtained with methods (a) 
and (e) of Table [2] normalised to the error bars, and compared 
with a Gaussian with the same area and cr = 0.33. 

While it is certainly possible to argue that certain algorithms 
for data reduction lead to more robust results than other meth- 
ods, some alternative choices may be equally reasonable. The 
fact that there are alternative and reasonable choices during the 
process of data reduction adds a contribution to the error bar of 
the mean longitudinal field. This contribution is difficult to quan- 
tify in general, but it is a non-negligible fraction of the error bar 
from photon noise. 

The results discussed in the remainder of this paper are based 
on a full reduction of the entire FORS 1 "fast mode" data archive, 
making the following choices: no fiat-fielding is implemented, 
spectra are extracted with average extraction using a 12 pixel ra- 
dius aperture, and rebinned at the dispersion value; the Py pro- 
files are rectified using a 150 pixel wide Fourier filter, clipped 
with k = 3 (including also the two bins that are adjacent to 
the deviant one) as explained in Sect. 13.21 the field is calculated 
from the full spectrum, although the first 3 % and the last 3 % of 
spectral points are ignored, and setting the Lande factor = 1 in 
proximity of the H Balmer lines, and 1.25 everywhere else. This 
method is referred to as method (a) in Table|2] The choices made 
are relatively conservative ones, tend to alter the actual data as 
little as possible, and attempt to describe the uncertainties real- 
istically. 

We note that no rectification is implemented in the FORS 
pipeline, since it would bring the aim of the software to a high- 
level of data manipulation which is in fact outside of the scope 
of a common user tool. Therefore, rather than using the final 
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Fig. 4. The distribution of the difference between (B z ) values 
obtained with two slightly different algorithms for data reduc- 
tion, normalised to the error bars, and compared to a Gaussian 
with identical area, centred about zero, and with cr — 0.33. 



products of the pipeline, we have always extracted the individual 
beams and treated them individually according to the algorithm 
described above. 

3.6. Alternative techniques for field measurements 

Traditionally, two main techniques are used in the literature to 
measure the mean longitudinal magnetic field. Measurements 
can be obtained from low resolution spectropolarimetric data us- 
ing the least-square technique used in this paper, or (more com- 
monly) from high-resolution spectropolarimetric data treated 
with the Least-Square Deconvolution (LSD) technique. In the 
latter case, the mean longitudinal field is obtained from the mea- 
surement of the first-order moment of LSD Stokes V about the 
line centre, and of the equivalent width of LS D Stokes / - 
both profiles normalised to the continuum intensity (Donat iet al.l 
U997tlKochukhov et al.ll20loh . 

One may experiment with variations of these two techniques. 
Equation (O could be applied to high-resolution spectra, or to 
their corresponding LSD profiles. LSD profiles themselves may 
be obtained even from low-resolution data, allowing one to mea- 
sure (B z ) either with Eq. (0), or with the moment technique. In 
practice, none of these alternative methods has been adopted so 
far, and it still has to be demonstrated that one can extract fully 
meaningful LSD profiles from low -resolution spectra. While dis- 
cussing alternative techniques for (B z ) measurements is beyond 
the scope of this paper, we report on our LSD experiments 
with three stars that were all observed with both grisms 600 B 
and 1200B: HD 358 (a non-magnetic star), HD 201601, and 
HD 94660 (both magnetic). Using the moment technique, no 
field was detected in the LSD profiles of HD 358, while the fields 
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Table 2. Longitudinal magnetic fields of three stars obtained 
with different algorithms. 





HD 94660 


HD 96441 


HD 171184 


Method 


(B-> (G) 


<S ; > (G) 


<S ; > (G) 


(a) 


-2506± 68 


-74+ 39 


125+41 


(b) 


-2622± 63 


-65+ 39 


139+41 


(c) 


-2549± 69 


-32+47 


150+45 


(d) 


-2512+ 68 


-71+39 


166+40 


(c) 


-2452+ 69 


-78+ 39 


127+ 42 


(f) 


-2496± 68 


-79+ 39 


148+40 


(g) 


-2508± 69 


-77+ 39 


127+41 


(h) 


-2453+71 


-74+ 40 


110+ 62 


(J) 


-2499+ 67 


-75+39 


117+40 


(k) 


-2811+ 101 


-168+47 


171+56 


(1) 


-2406+99 


123+77 


53+62 


(m) 


-3005+ 146 


-185+57 


126+71 



Keys to the methods, (a): as described at the end of Sect. 13.51 (b): 
same as (a), but with flat-fielding correction; (c): same as (a), but with 
optimal extraction; (d): same as (a), but using a 6 pixel aperture; (e): 
same as (a), but using a 18 pixel aperture; (f): same as (a), but without 
clipping; (g): same as (a), but without rectifying; (h): same as (a), with 
a 250 pixel wide filter; (j): same as (a), trimming only the first and 
the last 10 spectral points; (k): same as (a), but using Balmer lines 
only; (1): same as (a), but using metal lines, (m): same as (k) (i.e., using 
Balmer lines only), but applying a 2 pixel rebinning; 
For all cases, we have computed error bars using Eq. ( lilt . All observa- 
tions were obtained with grism 600B and a 0.5" slit width. Observing 
dates were 2003-02-08, 2004-04-17 and 2003-08-29 for HD 94660, 
HD 96441, and HD 171184, respectively. 



1.1 - 



1 » 



p 
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grism 600B (R-1600) 




grism 1200B (i?~2700 
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v (km s~ ! ) 
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Fig. 5. LSD profiles for HD 94660 observed with grism 600 B 
(left panel) and 1200 B (right panel). Red solid line: Stokes V/I c ; 
blue dotted line: N/I c ; black solid line Stokes I/I c . V/I c and N/I c 
are offset to +1 % for display purpose. 



measured in HD 201601 and HD 94660 were found to be ~ 25 % 
smaller (in absolute value) and more dependent on the adopted 
grism, than those obtained with the technique of Sect. 13.31 In all 
three cases, we found no field detection in the null profiles, al- 
though for HD 94660 observed with grism 1200 B, the null field 



value was = 2.4 cr (while with the least-square technique, the null 
field was consistent with within 1.6 cr). We performed numer- 
ical simulations, which confirmed that the (B z ) values obtained 
with the moment technique applied to LSD profiles tend to un- 
derestimat e the actual field s trength, at least for <: 1 kG field (see 
also iKochukhov et al.l 12010). Photon-noise error bars obtained 
with the two different methods were found similar among them- 
selves, but for the field obtained from LSD profiles we found an 
additional large uncertainty due to the continuum normalisation: 
the (B z ) values measured with the moment technique applied 
to the LSD profiles would change by more than 1 cr if spectra 
were normalised according to different but still reasonable crite- 
ria. (We note that normalisation criteria of single order spectra 
are somewhat more arbitrary than in case of echelle spectra.) 

In spite of these drawbacks, we note that in fast rotating stars, 
LSD V profiles may potentially reveal the signature of a mag- 
netic field even when its average longitudinal component is zero. 
Therefore, applying LSD techniques to FORS data may be still 
of some interest. Figure|5]shows the LSD profiles of HD 94660, 
observed with grisms 600 B and 1200B. 



4. Quality checks of the field measurements 

Inspection of the FORS 1 archive data shows that between 2000 
and 2009 about 1000 hours of telescope time were dedicated to 
acquiring stellar circular spectropolarimetric data for magnetic 
field measurements. Almost 600 hours of telescope time were 
granted in service mode, and 42 nights were granted in visitor 
mode. 

Figure [6] gives a qualitative technical overview of the na- 
ture of the observations, i.e., the distributions of the spectral 
resolution, of the total exposure time per observing series, of 
the peak signal-to-noise ratio per A, and the signal-to-noise ra- 
tio integrated over all pixel bins used for the field determi- 
nation (when using the full spectrum). The former two his- 
tograms demonstrate that the FORS instrument was often used 
for very bright objects, pushing the spectral resolution close 
to the limit fixed by sampling, and with substantial slit losses. 
This suggests that some FORS 1 users would have rather benefit- 
ted from high-resolution spectropolarimetric facilities at smaller 
telescopes. The latter two histograms show that most of the ob- 
servations were obtained with a ultra-high signal-to-noise ratio, 
hence FORS field measurements are generally close to the high- 
est possible precision obtainable with a reasonable use of the 
telescope time. The actual precision of the FORS field measure- 
ment is the subject of this Section. 

The discussion of Sect. [33] shows that using different reduc- 
tion procedures to infer magnetic field strengths may lead to sig- 
nificantly different results, in terms of both the mean longitudinal 
field (B z ) and its associated uncertainty cr^-) obtained from any 
specific observation. In this Section we fix the data reduction al- 
gorithm, and explore how the deduced fields and uncertainties 
still lead to a number of inconsistencies. This discussion has two 
goals: (1) to use what we know about some of the stars observed 
to try to learn more about the characteristics of the FORS 1 mea- 
surements, and (2) to examine the reliability of various results, 
particularly announced discoveries of new magnetic stars, based 
on FORS 1 measurements. 



4.1. Internal tests 

In principle, calibration and testing of the measurements from 
any spectropolarimeter should rely on calibration observations 
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Table 3. (B z ) measurements of a And = HD 358 




1000 2000 3000 4000 600 1200 1800 2400 3000 3600 4200 

Spectral resolution at central wavelength Exposure time (sec) 




2000 4000 6000 8000 10 s 2xl0 5 3xl0 5 

Peak signal-to-noise ratio per k Effective signal-to-noise ratio per A 

Fig. 6. Some characteristics of the FORS spectropolarimetric 
observations retrieved in the archive. 



made on a set of stars whose magnetic fields are well charac- 
terised by previous work, particularly including some which are 
known to be non-magnetic at the level of precision expected for 
the instrument. Unfortunately, such calibration measurements 
were not systematically obtained during the working life of 
FORS 1 . Therefore we have searched the catalogue of measure- 
ments for observations of stars with well known magnetic char- 
acteristics to test the actual precision that has been achieved with 
the FORS1 observations. 



4.1 .1 . Repeated observations of well-known magnetic stars 

A dispersion test that has been very useful in determining how 
realistic the measurement uncertainties are of other instruments 
such as the MuSiCoS spectropolarimeter (e.g. JWade et al.l2 000) 
is to examine the measurements of (B z ) as a function of the star's 
rotation phase around the mean variation. Thus, for example, 17 
measures of (B z ) of B CrB with MuSiCoS with uncertainties of 
20 to 30 G were fully consistent with the mean curve at signal 
to noise ratios of up to 30: 1 . This kind of instrument monitor- 
ing was not performed with the FORS 1 instrument, and we sug- 
gest it should be implemented in the FORS2 calibration plan. A 
somewhat useful test was performed by some users who repeat- 
edly observed the magnetic star HD 94660. The star's rotation 
period is not known, but is probably longer than two years. The 
star's longitudinal field appears nearly (but not perfectly) con- 
stant and about -2 kG, thus spectropolarimetric observations of 
HD 94660 allow one to assess if the polarimetric optics are at 
least approximately aligned. 



4.1 .2. Repeated observations of non-magnetic stars 

One particular test that should be carried out frequently as part 
of routine calibration of any spectropolarimeter used for mag- 
netic measurements is the observation of stars known to be non- 
magnetic, or at least to be non-magnetic to within rather small 



Date 


Hubris et al. (2006b) 


This work 


yyyy-mm-dd 


(B z ) (G) 


(B z ) (G) 


<N Z ) (G) 


2003-09-28 


-261+73 


-296+ 121 


-51+114 


2003-11-20 


12±82 


160± 104 


-143+ 121 


2005-05-29 


-109±49 


-40± 113 


216+ 103 


2005-09-16 


-73+20 


23+33 


12+31 


2005-09-17 


-30+ 30 


-12+39 


-30± 34 


2005-09-25 


-108±23 


-36± 37 


-26± 36 


2007-11-28 




107± 34 


1 1 1± 31 


2007-11-29 




-20± 24 


-29+ 20 



uncertainties (a few G or tens of G). Such measurements are a 
particularly powerful test of the uncertainties computed for indi- 
vidual observations, since the actual field value expected (zero) 
is known with high preci sion. An example o f such a test was the 
monitoring of 6 Eri by iDonati et al] (119971) with the SEMPOL 
instrument at the Anglo-Australian Telescope (AAT). No main 
sequence stars known with high precision to be non-magnetic 
have been deliberately observed with FORS1. However, some 
of the peculiar A stars of the HgMn type, and in particular 
a And = HD 358, have been repeatedly observed, probably be- 
cause they were suspected to be mag netic. Numerous investiga- 
tions of magnetism in such stars (e .g. Borra & Landstreet 1980; 
Glagolevskii et al.ll985tlShoriin et al.ll2002tlAuriere et al]2 010; 
Mak aganiuk et al.ll201 ll) have found no convincing evidence for 
fields in this class of stars with uncertainties often in the range 
of a few G up to a few tens of G. Furthermore, a And has 
been specifically investigated with the MuSiCoS (13 observa- 
tions, typical cr ~ 30 - 60 G) and ESPaDOnS (5 observations, 
typical cr ~ 6 — 19 G) spectropolarimeters without any signifi- 
cant detection of a non-zero longitudinal field or of any S tokes 
V signature indicative of a non-zero field (IWade et aUl2006l) . 

The field of a And has been measured eight times 
with FORS 1 . The fi rst six observations have been published 
dHubrig et al.ll2006bl) . Table [3] lists the published (B z ) measure- 
ments and the (B : ) values of the same observations resulting 
from our reduction. 

We note several interesting features of the data in this table. 
First, our new reduction invariably gives somewhat larger mea- 
surement uncertainties than the published data. Secondly, the 
published data report three field detections at the level of 3 cr or 
higher; our new reduction finds only a single detection at the 3<x 
level, and it is not one of the three initially reported 3 cr measure- 
ments. Furthermore, the null field (N z ) (i.e., the magnetic field 
calculated using the null profiles instead of the Py profiles) is 
also detected once at about the 3 cr level. Thus, our new results, 
independently of any external evidence, make the detection of a 
longitudinal field in this star from FORS 1 data very uncertain. 
The quite substantial differences between the field values pub- 
lished previously and those we have measured are probably not 
due to errors of reduction, but to the sensitivity of the results to 
the specific reduction method used to obtain the (B z ) value. 

If we now assume that t he more precise measurements re- 
ported by Wad e et alj d2006l) establish that no longitudinal field 
is present in a And at the level of 20 or 30 G, so that effectively 
we assume that a And is a null standard, we may regard the 
results of Table [3] as providing useful information about the be- 
haviour of the (B z ) values and associated uncertainties obtained 
from FORS1 measurements. In particular, both (B z ) and (N z ) 
values scatter from zero by more than is expected from the stan- 
dard deviation of their individual measurements. Our results sug- 
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gests that even the "external" error of Eq. (TTTt . which is still 
"internal" to a single measurement, leads to occasional underes- 
timates of the actual (B z ) error bars. 



4.1 .3. Time series 

A further test of the stability of measured field values, and of the 
associated uncertainties, can be obtained by studying the statis- 
tical properties of a long series of magnetic observations of the 
sa me star taken during a single night. One example is the study 
by Hubrig et al. (2004b), who observed six known cool magnetic 
rapidly oscillating Ap stars (roAp stars). These observations had 
a cadence of the order of one frame per minute, and were typi- 
cally continued for one to a few hours over ten or twenty pulsa- 
tion cycles. The observations were intended to search for possi- 
ble periodic variations in the magnetic field (B z ) at known pulsa- 
tion periods. According to the published result of this study, no 
short-period field variations were reliably detected. If we make 
the assumption that during the short time interval in which each 
star was observed, the field was in fact constant, we can compare 
the standard deviation of the values of (B z )j about their mean 



Table 4. Null field values and their error bars for six roAp stars 



Star 


ii 




cr ( N; ) 




(TV;)'"' 


w 


e 


HD 83368 


31 


-5 


44 


67 


-1 


8 


12 


HD 101065 


16 


-9 


40 


37 


-6 


11 


9 


HD 128898 


14 


11 


53 


101 


-2 


15 


27 


HD 137949 


16 


8 


41 


70 


17 


11 


17 


HD 201601 


17 


86 


80 


130 


72 


19 


32 


HD 217522 


50 


39 


110 


190 


23 


16 


27 


HD 101065 


36 


4 


45 


58 


-6 


8 


10 



Then we have calculated: 

the weighted mean of the (A^) (,) value 



(N z ) 



I?W> ffl /«>) a 



the average "internal" standard deviation 



value (B z ) with the computed uncertainties cr^ B , 
deviation (the "external error") is similar to the published val 
ues of crj'^ ^ ("internal errors"), then the comparision supports 
the correctness of the internal errors. Excess standard deviation 
indicates problems with deduced values of cr< B >. 

This comparison was car ried out by the team who did the 
observations (see Table 9 in iHubrig et al.l 2004b). Of the six 
stars reported, only two had external errors close to the inter- 
nal ones. The other four stars all showed excess scatter, up to an 
extreme external error of more than three times the internal er- 
rors for HD 201601. One interpretation of this excess, of course, 
would be that these roAp stars were actually displaying pulsa- 
tional field variations, but in only one case could a frequency be 
identified that is similar to one of the known pulsation frequen- 
cies of the roAp star, and a second observing ru n on this star 
(HP 1 01065) revealed no signal at that frequency (Hubri get alJ 
l2004bl) . Furthermore, HD 101065 is one of the two stars ob- 
served for which external and internal errors are very similar 
in magnitude; none of the stars with a significant discrepancy 
between internal and external errors showed significant varia- 
tions at known pulsation periods. Thus it seems far more likely 
that the difference between internal and external errors reveals 
something about FORS 1 and/or the reduction methods used for 
spectropolarimetric observations. 

We have repeated the dispersion analysis for these stars, us- 
ing our new reductions^ However, rather than using (B z ), for 
this kind of analysis we have preferred to use the mean longitu- 
dinal field (A^) from the null profiles. The reason is that Eq. (0 is 
strictly valid only in the weak-field regime, and for non blended 
spectral lines. These conditions are not met in the observed roAp 
stars, hence it is more likely that the null field (N z ) has a more 
linear response to the photon noise than the magnetic field (B z ). 
We have proceeded in the following way. From all M frames of 
a given observing series we have selected 4 x n frames, where n 
is the integer part of M/4, and calculated n (N z )® + cr*^ , values. 



— l v (,) 
If this standard cr{N z ) - - / ,tr <A g , 



the "external" standard deviation 

^2^ 1/2 



cr, 



( z;m> (0 - 



n 



1 



and error of the mean 



{ Z" «N Z )® - (N z )) 



2\ 



1/2 



cr 



<iV-> 



.1/2 



(12) 



(13) 



(14) 



(15) 



2 One particular feature of the newly reduced data is the change in 
sign of the data for HD 83368 compared to the earlier publication by 



Hubri g et al] {2004b). This point has been carefu lly investigated, and 
we believe that the signs in Hu brig et alJ i2004bh are incorrect for all 
the data for this star, and should be changed to positive signs. 



n(n - 1) 

We have also combined all An frames and obtained a (B z )^ ± 
cr'" ] , value for each star. 

Table|4]reports these quantities for each observed star. While 
there is good agreement between the (N z ) and (N Z ) M values, 
in all cases but one we found cr^ N ^ > crT^ and, consistently, 

e > cr^ . . Our data sets appear to reveal excess noise over what 
would be expected from the internal uncertainties. Excess noise 
arises irregularly: it increases the standard deviation by a fac- 
tor varying from 1 to about 2, qualitatively suggesting that the 
standard errors of measurement of (B z ) from internal dispersions 
within each measurement are underestimates of the true uncer- 
tainties. 



4.1 .4. The distribution of the magnetic field from the null 
profiles 

A final internal test is to compute the distribution over the entire 
FORS1 sample of field measurements of the null field values 
normalized by their error bars. 

The null profiles should represent measurements of zero 
fields, as the combinations of different waveplate settings used 
are chosen so that the real polarisation signal cancels out, to the 
precision of the measurements. If the field values measured from 
the null profiles have a distribution determined by photon noise 
only, this distribution should closely resemble to a gaussian dis- 
tribution with cr = 1. Figure [7] shows that there is a significant 
excess of (N z )/o-^_) points outside of the interval [-3,3]. The 
standard deviation of the (AQ/cr^) is 1.37, 1.20, 1.36 for the 
Balmer lines, metal lines, and the full spectrum, respectively. 
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Fig. 7. Histogram of the (N z ) values calculated from Balmer 
lines, and normalised to the error bars cr^ z ). The red solid lines 
show the gaussian with cr — 1, and the blue solid line a gaus- 
sian with <x = 1.365, both with the same areas as the observed 
distribution. The vertical lines mark the limits at -3 and 3. 



Fig. 8. The distribution of the differences between the null field 
values measured after flexure simulation and their original val- 
ues, normalised to the error bars, compared to a Gaussian curve 
centred about zero and with <x = 0.2. Null field values were cal- 
culated from the H Balmer lines. 



4.2. Possible explanations for the internal inconsistencies 

The results obtained for non magnetic stars, and the analysis of 
the time series and of the distribution of the null profiles, sug- 
gest that error bars, even when calculated with Eq. (fTTT) - which 
are generally higher than those obtained Eq. ( [Tol l, are underes- 
timated. Our conclusion is that Poisson noise is not the only 
source of uncertainty of the field measurements. We have already 
seen that an additional source of uncertainty is introduced by 
the choice of the algorithm for data reduction, but the results of 
Sect. 14.1.31 and I4TTT41 were obtained consistently using the same 
algorithm. Therefore we suggest that some physical, probably 
instrumental, limit is hit, especially when measurements are ob- 
tained with ultra-high signal-to-noise ratio. 

Possible reasons that could contribute to explaining this phe- 
nomenon are: 

i) Undetected problems with the raw data, such as saturated 
frames, or various kinds of errors in the organisation of the sci- 
entific frames and their calibrations, which range from the awk- 
ward mixing up of frames of different stars, to more difficult to 
detect hardware failures which could be responsible for an erro- 
neous recording of the position angle of the retarder waveplate 
in the fits-headers (these problems are known to have occurred 
at some time in the FORS1 data). 

ii) Change of the line profiles during the observations for reasons 
intrinsically due to the source (e.g., stellar pulsations). 

Hi) Changes of the radial velocity of the target occurring during 
the observing series. 

iv) Small instrument flexures occurring due to instrument motion 
during the observing series. 

v) For series of very short exposures in good seeing, changes in 
the mean position of the star in the slit from one exposure to 
the next, with corresponding changes in the mean wavelengths 



associated with each pixel. We note that many observations were 
obtained with exposure time of 0.25 sec per frame. 

We have clearly tried to minimise problems similar to those 
described in i), and issues ii) and Hi) can be investigated case by 
case. In this regard, it is interesting to note that many dubious 
field detections appear in slowly pulsating stars (see Sect. I5.4t . 

In the case of FORS 1 and FORS2, instrument flexures are 
kept at a minimum. However, we must note that ultra-high pre- 
cision polarimetric measurements can be performed only with 
an extremely stable instrument, and it is for that reason that the 
ESPaDOnS and NARVAL instruments are sitting in thermally 
and mechanically isolated benches. The question is what level 
of stability the Cassegrain mounted FORS instruments can guar- 
antee, and how possible instrument flexures would impact the 
(B z ) and (A^) measurements. This does not have a simple an- 
swer. Numerical simulations on a single spectral line show that 
a drift in wavelength equal in both beams (e.g., due to a star's 
radial velocity change) has a much stronger impact on the null 
profiles than on the Stokes profiles. A wavelength offset that acts 
differentially on the two beams, and that depends on the posi- 
tion of the retarder waveplate, has a larger impact on the Stokes 
profiles than on the null profilefl We note also that instrument 
flexures affect not only the Stokes and the null profiles, but also 
the normalised derivati ve of the flux (whic h is calculated as ex- 
plained in Sect. 4.2 of Bagn ulo et al.ll2002h . and, finally, that the 

" This case is illustrated in Bagn ulo et ail {2009), who also claimed 
that the null profiles obtained with the difference method are more sen- 
sitive to differential wavelength offsets than the null profiles obtained 
with the ratio method. In fact, we found that this claim is wrong. Both 
ways to calculate the null profiles are equivalent, and in the cases (c) 
and (d) presented in their Fig. 6, both null profiles obtained with the 
ratio and the different methods are zero. 
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impact of all these effects on (B z ) and (N z ) values depends also 
on line profile and blending. 

In an attempt to roughly evaluate the potential consequences 
of instrument flexures, we have performed the following test. 
In the data reduction process, we have numerically simulated a 
drift in wavelength continuously and homogeneously acting be- 
tween the first and the last exposure. The total amplitude of the 
drift (i.e., calculated between the last exposure and the first ex- 
posure), was equal to a quarter of the spectral bin. Then we have 
computed the difference A(N Z ) between the (N z ) value so ob- 
tained, and the value obtained without introducing any offset. We 
have repeated this simulation for all FORS archive datasets, and 
calculated the distribution of the ratio A(N z )/o"(ff-)- Figure[8]dis- 
plays the results obtained from H Balmer lines. This histogram 
shows that even a nearly negligible offset during the observations 
may produce changes in the (N z ) determinations comparable to 
the error bars. For the narrower metal lines, the consequences of 
flexures are more pronounced, and the (N z ) offsets being twice 
larger than when calculated from H Balmer lines. Qualitatively 
similar results was retrieved for the (B z ) values. However, we 
should note that for a certain dataset, a given simulated drift 
does not necessarily produce an offset with identical magnitude 
for both the (B z ) and (N z ) values. This suggests that, at least 
for small flexures, the warning coming from the (N z ) detections 
has a statistical relevance. For individual cases, a detection in 
(N z ) does not necessarily invalidate a (B z ) detection in the same 
dataset. Viceversa, a (B z ) detection may be spurious even if there 
is no detection in the null profile. 

According to the FORS instrument manual, in imaging mode 
the instrument is stable within 0.2 pixel during 1 hour exposure. 
However, Fig. 2.4 of that manual shows that flexures as large 
as a quarter of a pixel are not uncommon if the instrument ro- 
tator rotates by 10 - 20°. Furthermore, the instrument manual 
refers to a situation in imaging mode; the additional optical el- 
ements introduced in spectro-polarimetric mode (PMOS mask, 
slit mask, Wollaston prism, retarder waveplate, grism, and pos- 
sibly a filter) certainly cannot improve the instrument stability. 
While we conclude that instrument flexures may explain at least 
in part why FORS field measurement uncertainties seem to be 
underestimated, we must also point out that we were unable to 
find a correlation between (N z )/cT(n z ) values and the airmass. 



5. The incidence of magnetic fields in various 
classes of stars 

Our conclusion that the error bars computed from photon statis- 
tics underestimate by at least 30 %-40 % the real uncertainties 
of field measurements has significant implications for the body 
of data obtained from FORS1, in particular for those measure- 
ments which have been used to announce new discoveries of a 
magnetic field at a relatively low significance level. Because of 
problems that spurious discovery announcements can cause, in 
the rest of this paper we concentrate on those stars in which a 
magnetic field was detected with FORS 1 at $ 6 <x level, and 
present a detailed comparison of our new results with previously 
announced discoveries of magnetic fields. A homogeneous and 
comprehensive catalogue of all FORS1 measurements will be 
published in a forthcoming paper. 

The results of our comparison are summarised in Table [5] 
The stars are grouped into families of various kinds, identified 
by bold-face titles such as "HgMn stars", "Be stars", etc. The 
first two columns in this table list the star name, and the mean 
MJD of the observation as computed for the frames used in our 



data reductions (a blank in one of these two columns means that 
the line refers to a second published value of the same data as ap- 
pears immediately above). The third and fourth columns give the 
value of (B z ) and its uncertainty as published, and a code denot- 
ing the reference where the field measurement was published. 
The fifth column reports the (B z ) value and its uncertainty ob- 
tained from our adopted reduction process. The sixth column 
includes a code character describing our conclusion on whether 
a field is detected in the star: D means that the star is definitfy 
magnetic; P means that a field is possibly present, but that the 
currently available measurements do not provide firm evidence 
for that; N means that none of the measurements shows evidence 
of the presence of a magnetic field. Note that a few of the entries 
in Table [5] have MJD values that differ from those published, 
mostly by 0.5 day. Our values are computed directly and uni- 
formly from the data in the fits headers of the observations, and 
should replace the published values where these differ. 

Note that Table [5] does not include the "marginal" fields 
det ections in Ap and Bp stars - most of which are reported 
byjHubrig et al. (2004b), Bagnul o et all d2006l). iHubrig et all 
(2006bJ) and lKochukhov & Bagnulo! (l2006h . It is already clearly 
demonstrated that magnetic fields are detected in essentially ev- 
ery A p star if high enough precision is reached dAuriere et alj 
2007). Therefore, the reliability of the detection of weak fields 
reported in these papers will not greatly change what we know 
about the occurrence of magnetic Ap and Bp stars. Accordingly, 
we postpone a discussion of these data to our forthcoming paper. 

The uncertainty of our field measurements of Table [5]is cal- 
culated with Eq. (fTTT i. which is still an underestimate of the real 
error bar. Unfortunately, the source of errors that are not due 
to photon noise cannot be estimated for the individual cases. 
Although one may be strongly tempted to multiply all error es- 
timates by a factor ~ 1.5, the correct way to proceed is to re- 
member that 3 cr level is not a safe detection threshold. For that 
reason, the discussion presented in this Section will also rely 
on checking whether a field detection has been repeatedly con- 
firmed with numerous FORS1 measurements, and on corrobo- 
rating or conflicting field measurements made with other spec- 
tropolarimeters. 

We found that several important FORS1 discoveries an- 
nounced in the past are simply not confirmed by our new data 
reduction, even without taking into account a possible underesti- 
mate of our error bar. In many cases, in contrast to the originally 
published result, the absolute value of our new field measure- 
ment does not exceed three times its associated cr^ B value esti- 
mated with Eq. (fTTI) . We have two possible general explanations 
for this finding. The first one is based on the general consider- 
ations of Sect. [33] that different data reduction methods lead to 
slightly different results. Also, in some cases we have removed 
some frames (e.g., because of saturation), or used a slightly dif- 
ferent frame selection than used for the original publication. In 
an attempt to keep this aspect under control, we have compiled a 
very large catalogue that includes, for each observation, all (B z ) 
values obtained with several different data reduction methods 
(e.g, without rectifying Stokes profiles, or using optimal extrac- 
tion instead of average extraction). In the following discussion 
we will explicitly comment on those cases in which one or more 
methods for data reduction alternative to the one used to compile 
Table|5]would lead to a qualitatively different result. These cases 
are comparatively rare. A second explanation of the discrepan- 
cies between our results and those published in the literature is 
that our error bars are generally larger than those obtained in pre- 
vious works. In the following we comment in detail on our new 
results. 
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5.1. Am, HgMn and PGa stars 

The HgMn peculiar late B stars have already been mentioned in 
Sect. I4.1.21 where we concluded that the overwhelming bulk of 
evidence is against the presence of fields having (B z ) values of 
more than a couple of tens of G. Table [5] contains six measures 
claimed as field detections of four such stars by Hubri get al.l 
(2006b). The new reductions have increased the standard errors 
of all the measurements by factors of 1 .3 or more, in two cases 
by much more. None of the measurements now shows a field at 
the 3 cr level, and three of the five measures that were reported 
to be over 100 G and detected at even 4 cr or more are now only 
some tens of G, and at roughly the 1 cr level. In contrast, two of 
the large values are still relatively large, but are no longer signif- 
icant at the 3 cr level, and may represent the occasional spurious 
large field values that appear in a few percent of FORS 1 field 
measurements. 

We conclude that the FORS1 data do not alter the general 
non-detection of magnetic fields in HgMn stars. There is still no 
significant evidence of fields in such stars. 

The PGa stars are considered to be the extension of the 
HgMn phenomenon to higher effect ive temperatures. A survey 
of several such stars some years ago dBorra et al.ll 1983b revealed 
no fields, with unce rtainties of hundreds of G. The large FORS 1 
survey reported bv IHubrig et al.l (l2006bl) included two PGa stars 
observed with much smaller uncertainties. One of these stars, 
HD 19400, was reported to have a field. Our re-reduction of 
these data show no significant field, and we consider this de- 
tection to be spurious. 

The Am (metallic-line) stars are considered to be analogues 
of the HgMn stars at lower r e ff. There have been only a few 
observations of such stars with FORS1, and no claims of field 
detections based on these observations have been made. 



5.2. Classical Be stars 

By classical Be stars we mean those B stars which show ev- 
idence (at least some of the time) of centrifugally supported, 
slowly evolving, fairly massive disks. For many years some the- 
oreticians have argued that magnetic fields may play an impor- 
tant role in the creation and/or maintenance of the circumstellar 
disks. However, there have not been very m any serious searches 
for magnetic fields in such stars. A search by Bar ker et"a"D (l985) 
with uncertainties in the range of 50 - 250 G did not reveal any 
fields. 

On the basis of data from the MuSiCoS spe ctropolarime- 
ter.lHenrichs et a il d2000h . INeiner et alJ d200l and lDonati et a . 
j200lF have reported the detection of a weak field in HD 205021 
= BCep, a double system composed of a Be star and a pul- 
sating star - the latter be ing the proto-type of BCep pulsating 



stars. ISchnerr et a fl (120061) have shown that the primary star of 
the B Cep system is the pulsating and magnetic star, while the 
Be star is a non magnetic companion that can be separated with 
speckle interferometry. T hus the detected field is not connected 
with the Be phenomenon. INeiner et alJ (l2003al) reported a mag- 
netic field in the classical Be star u> Oh, but it is not clear that 
the data provide a co nvincing case for the presence of a field. 
ISilvester et afl d2009l) have observed the two classical Be stars 
HD 148184 and HD 181615 several times each with ESPaDOnS, 
with uncertainties of order 100 G and 10 G respectively, but 
found neither significantly non-zero values of (B z ), nor signfi- 
cant polarisation signatures. Thus from the limited data available 
from instruments other than FORS 1 , there is no strong evidence 
for magnetic fields in any classical Be star. 



The largest recent surveys for fields in classical Be stars 
have been carried out using FORS1 bv IHubrig etail d20 07b) 
and IHubrig et al.l d2009bl) (summarised by lYudin et all 1201 ll), 
both s earching for magnetism in field Be stars, and by McSwain 
(2008), who studied normal, chemically peculiar, and emission- 
line B stars in the young open cluster NGC 3766. Between them, 
these surveys have included more than 40 Be and possible Be 
stars, of which 10 stars were reporte d to have weak but s ignif- 
icant magnetic fields . The reports by IHubrig et alJ d2007bl) and 
IHubrig et all (l2009bl) suggest that weak (~ 100 G) fields may be 
ubiquitous in Be stars, a result that would have important conse- 
quences for modelling the Be phenomenon. 

Note that there is some doubt as whether all of the stars in 
these surveys are Be stars. For example, NGC 3766 MG 94, 111, 
and 17 6 are listed as a H e strong star, normal star, and normal 
star by McSwain (2008 ), but a s B stars with emission lines in 
Table 2 of Hubrig et al. (2009b). In Table[5]we followed the clas- 
sication bv iMcSwainl d2008l) . Furthermore, there are differences 
among these pu blications as to which stars show field detections; 
IMcSwainl d2008h did not find a significant field in the normal B 
stars NGC 3766 MG 1 1 1 or NGC 3766 MG 176, nor in the Be 
star NGC 3766 MG 200, all o f which are reported , using the 
same data, as new detections by Hubrig et"aT]d2009bl) . 

When we compare the new reductions of these observations 
with the published ones, we see that only one of the reported 
(B z ) values that differ from zero by more than 3 cr are still signif- 
icantly non-zero. In part this is because most (but not all) of the 
measurement uncertainties resulting from our re-reductions are 
larger than those originally reported, and in part because most 
of the new (B z ) values are also smaller (in absolute value) than 
originally reported. The only remaining apparently significant 
detection is in HD 181615; however, the photon-noise error of 
the measurement that still shows field detection in our reduction 
is certainly very small, and the presence of a field of the claimed 
size is clearly contra dicted by the three negative observations of 
ISilvester et all d2009h . 

There is a huge discrepancy between the uncertainty of the 
published field detections for HD 148184 and our uncertain- 
ties in Table [5] These spectra are heavily polluted with emis- 
sion lines, and we were una ble to repro d uce the very tiny field 
uncertainties published by Hubri g~et alJ (12007b). However, we 
must note that the quality of our reduced data corresponding to 
the second detection is particularly poor, probably due to image 
quality. We note that the publis hed detections are co ntradicted 
by four negative observations of Silve ster et al.l (120091) with un- 
certainties mostly of about 70 G. 

Our reductions using the ESO pipeline give apparently sig- 
nificant field detections for an observation of HD 56014 and two 
observations of HD 209409, in which the original observers did 
not find report any significant fields (for that reason, the corre- 
sponding entries are not included in Table|5]). We also obtained a 
15 cr detection (!) in HD 224686, in which the original observers 
measured (B z ) = 74 + 24 G. All these data were obtained in the 
same two nights, MJD=54432 and 54433. We are very scepti- 
cal of our reductions, since they all show a significant signal in 
the bluest Balmer lines both in the Py and in the Ny profiles, 
which we believe the symptom of a not yet understood prob- 
lem with calibrations or with our data reduction. In Table [5] for 
HD 224686, instead of our field estimate, we left an empty space. 

Our global conclusi on is that most if not a ll of t he detections 
reported in Be stars by Hub rig et alJ d2007bl) and Hubri g et al. 
(2009b), as well as the marginal detection reported by Mc Swain 
d2008l) . are probably spurious, and that magnetic fields much 
above 100 G rarely if ever occur in classical Be stars. 



13 



S. Bagnulo et al.: Magnetic field measurements and their uncertainties: the FORS1 legacy 



5.3. Herbig AeBe stars 

The Herbig AeBe stars are pre-main sequence stars which are 
in the final stages of contraction, and which are destined to be- 
come main sequence A- or B-type stars. These stars are identi- 
fied by their locations in regions of current or very recent star 
formation, by the presence of optical emission lines (especially 
Ha) produced by circumstellar or accreting material, and by the 
presence of an infrared excess due to cool circumstellar material 
not yet accreted onto the star. The stars themselves have spectral 
types A or B. 

The first discov ery of a fiel d in a H erbig Ae star, HD 104237, 
was reported by iDonati et al.l d 1997b from observations at the 
AAT with an experimental spectropolarimeter. The star was ob- 
served again with SEMPOL at the AAT by E. Alecian (pri- 
vate communicatio n) who confirms the magnetic detection pre- 
viously reported by Dona ti et all (1 19971) . 

Based on FORS1 observations, the discovery of a mag- 
netic field in the He rbig Ae star HD 139614 was reported by 
Hubrig et all d2004al) . who measured (B z ) = -450 + 93 G. 
Wade et at] d2005l) re -analysed the FORS 1 dataset acquired by 
Hubrig et al.l d2004al) . a nd obtained the m ore marginal result 
(B z ) = -150 ± 50 G. IWade et all d2005l) reported also that 
ESPaDOnS data did not show any presence of a magnetic field. 
A second detection in HP 1396 14 «£,) = -116 + 34 G) was 
claimed bv lHubrig et al.l d2006d) . who also confirm e d their pre- 
vious measurement of -450 G. Later, iHubrig et al.l d2007al) re- 
vised both FORS 1 observations, and reported (B z ) - 1 12 ± 36 G 
and (B z ) = -93 + 14 G, for the first and the second measure- 
ment, respectively. A magnetic fi eld was not detected in a third 
FORS1 measurement reported by Hubrig et al. (2009c). Our re- 
reduction of all three observations clearly shows that no field is 
detected in any of these three measurements. This star has also 
been observed three tim es using the high-resolut i on spectropo- 
larim eter ESPaDOnS (IWade et alj|2005l l2007bt lAlecian et ail 
120111) . ESPaDOnS data, with (B z ) uncertainties of about 12 G, 
show no trace of magnetic field, and demonstrate quite definitely 
that no field anywhere near the level originally reported is de- 
tected in this star. Our conclusion is that the reported discovery 
of a magnetic field in HD 139614 is spurious. 

The results of a survey with FO RS1 of abou t 50 H erbig 
AeBe stars were briefly described by Wa de et ail d2005l) . The 
same paper also reported first results from the use of the new 
ESPaDOnS spectropolarimeter at the CFHT to search for fields 
in Herbig AeBe stars. From the FORS1 survey, magnetic fields 
were reported in two stars: HD 72106A (a very young main se- 
quence star whose companion, HD 72106B, is a Herbig star), 
and HD 101412. Although the uncertainties of the measurements 
of these two stars have increased in our re-reduction of these 
data, both fields are apparently still detected. Since detections 
are only at a significance level between 3 and 4 cr, FORS 1 data 
by themselves do not stron g ly estab lish the presence of magnetic 
fields. In fact, Wa de et all d2007al) revisited their FORS1 mea- 
surement of HD 72 106 A, finding that their data could not sup- 
port a field detection. The field of HD 72 1 6A was confirm ed by 
the ESPaDOnS component of this survey dWade et al.l2005l). and 
the sta r has subsequently been studied in detail bv lFolsom et alJ 
(2008|). T he presence of a field in HD 10 1412 has been fully con- 
firmed bv lAlecian etal] d2008allR l2009ah with SEMPOL obser- 
vations, and bv lHubrig et alJ d2009d) with further FORS1 field 
measurements at a level ^ 6 cr (therefore these measurements do 
not appear in Table [5]). T he field of HD 101412 has been studied 
in more detail by Hubrig et al.l d201 lal) . 



In the survey reported by IWade et all d2005l) . the star 
V380 Ori was observed with FORS 1 without a significant detec- 
tion. However, the same paper reports a clear detection of this 
star (on a different date) with ESPaDOnS. The magnetic field 
has been detected at more th an the 3 cr level a nu mber of times, 
and was studied in detail bv lAlecian etal] d2009bl) . 

Further s urveys of Herbig AeB e sta rs using FORS1 wer e 
described by IHubrig et all d2006d) and IHubrig et all d2007iJ) . 
Altogether, these papers reported the discovery of a mag- 
netic field in three Herbig stars: HD 31648, HD 144432, and 
HD 144668. The obse rvations of HP 3164 8, presented as a 
ne w field discovery by Hubrig et alJ (l2006d) . were re-reduced 
by IHubrig et al.l d2007al) . who reported a final (B z ) measure- 
ment below the 3 cr significance. For HD 144432, both pa- 
pers rep orted a similar field value, but the uncertainty pub- 
lished by Hubrig et all(l2007al) was a factor of almost 2.5 smaller 
than estimated by Hubri g et al.l (l2006d) . Two m easurements of 
HD 144668 published by Hubrig et al. (2006c) are m arginally 
significant, but do not exceed 3 cr. IHubrig et al.l (l2007al) revised 
both measurements, and found for this star a new field detec- 
tion with a 4cr sig nificance. We also note that a survey by 
IHubrig et all d2009d) . that we will review later in more detail, 
reports for HD 144668 a second, marginally significant detec- 
tion. Our re-reduction of all thes e observations bring s the field 
measurement of HD 31648 bv IHubrig et al.l (l2006c[). and the 
field measurement of HD 144668 by Hubrig et al.1 (I2007ab to 
zero wi thin errors. Our re-reduction confirms the field detec- 
tion by IHubrig et al.l d2006cl l2007al) in HD 144432 at almost 
5 cr level, as well as the ~ 3 cr detection f rom th e observations 
of HD 144668 obtained by IHubrig et al] d2009d) . These same 
three stars were also obs erved with FO RS 1 (with simil a r uncer - 
tainties to those of IHubrig et al.ll2006cl) bv lWade etai] d2007d) . 
who did not detect any significant fields. In addition, magnetic 
field measurements of HD 31648, HD 14 4432, and HP 14466 8 
have been obtained with ESPaPOnS by lAlecian et all d2011l) . 
These new measurements have error bars between 40 G and 
145 G, and do not show any strong indication of the presence 
of fields (although one field measurement for each star reaches 
2 cr). Thus, it appears that the reported detections of HP 144668 
and HP 144432 might possibly be real, but these detections are 
quite uncertain, and need substantial new evidence to be con- 
vincing. 



The FORS1 survey briefly described by IWade et al.1 (120051) 
was analysed in more detail bv lWade et al.l d2007a ). Altogether, 
68 field measurements of 50 stars were carried out. The data 
were re-reduced by the authors, who concluded that HP 101412 
and BF Ori probably host magnetic fields, while HP 361 12 and 
CPP -53 295 might be magnetic, but that these marginal detec- 
tions would require confirmation. After our new reduction, mag- 
netic fields are still detected in HP 101412 and CPP -53 295 
at a slightly more than 3 cr confidence, while the measurements 
of HP 361 12 and BF Ori no longer show field detection. As dis- 
cussed above, the detection of HP 101412 was later fully con- 
firmed. A second m easurement of CPP -53 295 (published by 
Hubrig et al. 2009c) was reported as field detection but does not 
appear significant in our re-reduction of the data. The only fur- 
ther data available for HP 361 12 or BF Ori are two observations 
with ESPaPOnS or NARVAL for eac h of the stars, which show 
no hint of fields (lAlecian et al. I I20TT1) . We conclude that a field 
might be present in CPP -53 295, but this result would cer- 
tainly need further confirmation. The marginal fiel d detections 
of HP 36112 and BF Ori reported bv IWade et ail d2007al) are 
probably spurious. 
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The mos t recent FORS1 survey of Herbig stars is due to 
IHubrig et al.1 (l2009d) . Nine field detections at more than the 3 cr 
level are reported, including two field detections in HD 101412 
and one in HD 190073, two stars previously known to be mag- 
netic Herbig AeBe stars. New field discoveries are claimed for 
six stars (on the basis of one measurement each), although one 
of these new magnetic stars is CPD-53 295, which was al - 
ready reported to be possibly magnetic by Wad eltai] d2007al) . 
Four of these six detections, including the measurement of 
CPD-53 295, disappear in our re-reduction, while two of them 
remain signficant. As discussed above, a 3 cr detection of a 
field in HD 144668, might be real but requires confirmation. 
The detection at the 5 - 6cr level of a field in HD 150193 is 
even stronger after our re-reduction than originally published. 
However, a single ESPaDOnS obs ervation (with an uncertainty 
of about 120 G) detects no field dAlecian et alJl201lh . 

Finall y, we note that our new reduction of the datasets ob- 
tained by Wad e et al.l (|2007a) yields two field detections in the 
star HD 97048 (which were not fo und in the original re ductions), 
at just slightly over the 3 cr level. IHubrig et al. (l2009d) also has 



one detection at about the 4 cr level, which by constrast is not 
retrieved in our new reduction. It is possible that this star has a 
weak field of a few hundred G, but we do not consider the evi- 
dence for this to be very strong. 

In summary, FORS1 observations have successfully been 
used to discover the field of one Herbig star, HD 101412, and 
of the binary com panion of another Herbig star, HD 72660A 
dWade et aLlBOOl , FORS 1 observations have shown that fields 
are possibly present, but certai nly not yet definitely established, 
in HD 97048 and HP 14466 8 dHubrig et al.ll2006cl). HD 144332 
(IHubrig et alj|2006d l2007al). CPD -5 3 295 dWade et alj|2007 at 
IHubrig et al.ll2009d) and HD 150193 dHubrig et aUl2009d) . All 
other published FORS 1 field discoveries of Herbig AeBe stars 
appear to be spurious. Detectable magnetic fields appear to oc- 
cur in only a few percent of all Herbig AeBe stars observed 
with FORS1. This result is consistent with the 7 % found by 
Aleci an et al.l d201 lb from the comparable size ESPaDOnS sur- 
vey. 



5.4. B Cephei pulsators and Slowly Pulsating B stars 

The /3 Cep pulsators are early B stars, with spectral types mostly 
earlier than B3, and masses in the range of about 10 to 20 M©, 
which pulsate in at least one radial p mode with periods in the 
range of 0.1 to 0.5 d. The Slowly Pulsating B stars (SPBs) are 
stars which are B2 or later, have slightly lower masses than (3 Cep 
stars, and show non-radial g mode pulsations with periods of or- 
der 1 d. It appears that these pulsations are excited by the kappa 
mechanism operating in a deep region of high iron opacity. 

A field was discov ered in the prototyp e B Cep star, B Cep 
= HD 205021 itse lf dHenrichs et all 120001: INeiner et alJ l200lt 
iDona ti et al. 2001). The field of this star appears to be very se- 
curely detected (see Sect. 15.2b . 

A magnetic field was reported in the B Cep star V2052 Oph = 
HD 163472 bv INeiner et al] d2003bh . again based on MuSiCoS 
observations. The measurements presented appear to support the 
claimed discovery, but are very close to the limit of detection. 
The reality o f this field is now clearly confirmed by NARVAL 
ob servations ( NeineretaLll2011l) . 

INeiner et aL ( 2003d) reported the presence of a magnetic 
field in the SPB star £ Cas = HD 3360 (still based on MuSiCoS 
data). Although these observations appear consistent with the 
presence of a field, they do not support fully convincing evi- 
dence for it. However, the presence of a magnetic field in ( Cas is 



definitely confirmed by further (unpublished) observations with 
MuSiCoS and NARVAL (C. Neiner, private communication). 
Major surv eys of both B Cep an d SP B pulsators have bee n 



Major surveys ot botn p Cep ana t>rti pulsators nave been 
carried out by iHubrig et all (l2006al) and IHubrig et al.l d2009al) . 
Nearly 70 stars were observed for magnetic fields using the 
FORS1 spectropolarimeter (19 known and suspected B Cep 
stars, 50 known and suspected SPBs). In this survey, fields were 
reported to have been detected in five B Cep stars and 26 SPB 
stars. If these results were confirmed, they would be quite im- 
portant; such high incidence of detected fields would strongly 
suggest that magnetism is intrinsically connected with the pulsa- 
tion phenomenon in early B-type stars, as it is in the cool rapidly 

osc illating Ap (roAp) sta rs. 

IHubrig et ail d2006al) and lHubrig et al.l d2009al) have reported 
discovery of a field in the B Cep star CMa = HD 46328, 
based on altogether 13 detections at approximately the 6 or 7cr 
level. The field always appears to be close to +350 G. These 
detections are furthe r supported by a single field detection by 
ISilvester et all d2009l) using ESPaDOnS at about th e 30 cr level 
by tw o ESPaDOnS field measurements reported bv lShultz et al.l 
(1201 lh . and by further ESPaDOnS observa tions not yet pub- 
lished in detail dFourtune-Ravard et alJl201 ll) . The re-reductions 
of the FORS1 data (not shown in Table [5] because the measure- 
ments exceed the 6cr upper limit for inclusion) also confirm 
the presence of a field. This discovery appears very robust, and 
clearly indicates the capability of FORS 1 to detect rather modest 
magnetic fields in B-type stars. 

The reported field detections in B Cep stars in Table [5] have 
mostly become insignificant in the new reductions. In particu- 
lar, the six reported detections of a field in HD 16582 all have 
decreased below the 3 cr significance limit, althou gh one mea- 
sureme nt not originally claimed as detection by Hubri eet al.l 
d2009al) . on MJD=54343.259, has risen to the apparently sig- 
nificant value of -104 + 19 G. A single ESPaDOnS field mea- 
surem ent detects no trace of a field with (B z ) uncertainty of 
10 G dSilvester et al.ll2009l) . While it is still possible that a field 
might be present in HD 16582, this must certainly be regarded 
at present as - at best - extremely uncertain. Only two pub- 
lished field detections, one each for HD 74575 and HD 136504, 
are still significant at the 3 cr level, and each of these measure- 
ments is only barely significant. Three field measurements with 
ESPaDOnS of HP 74575 show no trace of a field at the 6 G level 
dShultz etal.ll201ll) : we regard these measurements as clear ev- 
idence that no field has been detecte d in this star . On th e other 
hand, a field is definitely detected bv lShultz et al.l d201 lh in two 
out of two measurements of HP 136504 (=eLupi ); these mea- 
surem ents clearly confirm the field discovered by Hubri g et al.l 
(2009a). All the remaining reported detections in B Cep stars ap- 
pear to be the result essentially of underestimated measurement 
uncertainties; currently there is no significant evidence for fields 
based on these measurements. 

Of the more than 40 fi eld detections reporte d in SPB stars by 
IHubrig et all d2006al) and lHubrig et al.l d2009ah . all but five have 
decreased to non-detections in the new reductions. 

For HP 281 14, the re-reduced field strength is considerably 
larger than the published value (and of the opposite sign), but 
significant at the 4 cr level. Field detection is obtained only on 
the higher order Balmer lines, which suggests that our detection 
is spurious. A single ESPaPOnS observati on finds no trace of a 
field at the 30 G level dSilvester et al.ll2009l) : it is presently not at 
all clear that a field has been detected in HP 281 14. 

For HP 74195, the published detection is no longer signifi- 
cant, but a field detection on MJP=54143.072, contradicted by 
the non-significant published field value, is now found at slightly 
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more than 3 cr significance, so that a field might possibly be 
present. 

In the new reductions, HD 53921 still has one significant 
field detection, although with the opposite sign to the published 
value; the other two published detections are no longer signifi- 
cant. The star is a visual binary with 1.4"separation; in less than 
optimal seeing conditions the two components cannot be easily 
disentangled (e.g., on MJD=53630.401 and 53631.408 the two 
spectra were not distinguishable.) Note that a field has also been 
detected in this star in two HARPSpol measurements recently 
obtained by one of us (O. Kochukhov), so the detection of a field 
in this star appears to be secure. 

The new reductions of HD 15251 1 still have two apparently 
very significant detections (but we note that for reasons probably 
related to bad image quality, our reductions of data obtained on 
MJD=54609 lead to poor results). 

For 16 Peg = HD 208057, the new redu c tions c onfirm the 
original detection published by Hubrig et al. (2006a), and also 
indicate a significant detection of -195 + 60 G at slightly over 
the 3 cr significance level for the other measurement of this star 
(for which the original published field was not a d etection). Note 
that a field has also been marginally detected by ISilvester et alj 
(120091) in a single ESPaDOnS measurement of HD 208057; the 
reality of this field has also be en much more extensively con- 
firmed bv lHenrichs etal] d2009l) . 

In % 1 % of the observations, the FORS pipeline failed to 
reduce the scientific frames. The dataset of HD 161783 obtained 
on MJD=53520.308 is one of these rare cases, and our entry in 
Table|5]is an empty space. 

We finally note that the field measurement o f HD 215573 
obtaine d on MJD=53193.321 and published bv IHubrig eTaT] 
(2006a) is obtained from frames that seem to have reached 
saturation, even though the ADU count is just around 50000. 
We found that the CCD gain was set such that the electron 
level would reach the full well capacity before ADC saturation. 
Accordingly, we have discarded this dataset, and left an empty 
space in the entry of Table [5] The same problem affects many 
observations of programme ID 073.D-0466, dedicated to the ob- 
servations of SPB and Bp stars. 

Based on the re-reductions of the FORS1 data, as well as 
other published material, we conclude that the securely detected 
fields at present among the p Cep pulsators are those of ft Cep, 
V2052 Oph, CMa, and HD 136504 (=e Lupi). Securely de- 
tected fields occur in the SPB stars £ Cas = HD 3360, HD 53921, 
and 16 Peg = HD 208057. p Cep star HD 16582, and SPB 
stars HD 74195 and HD 15251 1 may possibly have weak fields, 
but these apparent detections still need to be firmly confirmed 
or rejected. The great m ajo rity of field d i scoveri es announced 
by IHubrig et alj d2006al) and IHubrig et al.l d2009al) are spurious. 
Based on this conclusion, it appears that, as in other kinds of 
upper main sequence st ars, magnetic fields are relatively rare. 
Gru nhut& Wadd d201lh estimate that fields are present in 16 % 
of pulsating B stars, only marginally more than the incidence 
of magnetic fields among A stars. More importantly, there is 
currently no significant case for considering that the pulsation 
properties of these B pulsators are intrinsically connected to the 
presence of weak magnetic fields. 

5.5. 6 Scuti pulsating late A-type stars 

Delta Scuti stars are late A- or early F-type stars, on or near 
the main sequence, that pulsate by the kappa mechanism with 
periods of a few hours. The pulsations are a mixture of radial 
and non-radial modes. 



Although 6 Set stars occur in stars of the same mass and 
evolution state as cool magnetic Ap stars, there are no confirmed 
magnetic Ap stars that show 6 Set pulsations. From photometry 
with Kepler mission, several possible magnetic Ap 5 Set stars 
have been reported, although the spectral classification of these 
stars is still quite uncertain. One definit e case of a field in a pul- 
sation y Dor star has been discovered (Ba lona et al.ll201 ll) . On 
the other hand, there does not appear to have been any substan- 
tial survey of known 5 Set stars for magnetic fields, although oc- 
casionally surveys of bright A- and F-type stars hav e included 
one or two 6 Set stars such as Altair and p Cas dLandstreell 
1982; Monin et a 1. 2002). No field detections have emerged from 
these observations, although it would be interesting to carry out 
a larger survey of 6 Set stars. 

Using FORS1. iKurtz et all d2008l) have observed the 6 Set 
star HD 21190, f or which an Ap classification was proposed 
dKoen et al.ll2001l) . and have reported an apparently significant 
field of (B z ) = 47 ± 13 G, although they comment that this de- 
tection requires confirmation. Our re-reduction of this measure- 
ment increases the uncertainty only a little, but the result is that 
there is no significant detection. We regard this detection as spu- 
rious. We also need to add that a few Ap classifications of known 
5 Scuti stars are inconclusive, since none was supported by a de- 
tailed spectroscopic study. A quick inspection of archival UVES 
spectra allows us to conclude in particular that HD 21 190 is not 
an Ap star. 



5.6. Non-peculiar B-type main sequence and giant stars 

Magnetic fields have begun to be detected in B-type stars that 
do not show the obvious chemical peculiarities that signal mag- 
netic Bp or He-peculiar stars. The stars p Cep, £ Cas, and 
CMa have been discussed above in the context of pulsating 
B stars. Fields have also been reported in r Sco = HD 149438 
dDonati et all l2006bl) HP 36982 = Par 1772 and NU Ori= 
HD 37061 dPetit et alJl2008h . 

FORS 1 surveys i n cluding such stars have been reported by 
IHubrig et all d2006bl). IMcSwainl d2008l) . iHubriget all d2009ah 
and IHubrig et al.l ([2009b). As noted in Table [5] detections of 
fields have been reported in several stars. In our new reduc- 
tions, the only remaining measurement that shows a significant 
field is one of two observations of HD 5208 9, whose field ap- 
pears to be significant at about the 5 cr level. IMcSwainl (2008) 
reported a marginal detection in NGC 3766 MG 45, appar- 
ently contradicted by our estimate of Table [5] obtained from 
the full spectrum. However, from H Balmer lines, we measure 
(B z ) = -258 + 77 G. In conclusion, a field might be present 
both these stars, but would certainly require confirmation. The 
other reported detections resulting from these surveys are prob- 
ably spurious. 

In addition, Hubrig et al.l d2008al) report a number of mag- 
netic field measurements of the BOVp star 9 Car = HD 93030, the 
brightest star in the open cluster IC 2602. Five of these measure- 
ments are significant at ^ 3 cr level, as listed in Table|5j although 
IHubrig et al.l d2008al) do not conclude that the star is magnetic. 
Because the 26 measurements were all obtained within about 
80 min, in our re-reduction of the field measurements we have 
combined all measurements into a single average value. The av- 
erage measurement does not show any signficant field, and we 
conclude that FORS data do not reveal presence of a magnetic 
field. 
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5. 7. O-type stars 

O-type stars are relatively difficult targets for magnetic field 
measurements, because they have few strong spectral lines, high 
projected rotational velocities (~ 100 kms~'), and often fairly 
strong emission filling in some of the photospheric absorption 
lines. The first detections came when a field wa s detected in 
l Or i C, the brightest star in the Trapez ium, bv | Ponati et all 
(120021) . and in the Of ?p star HD 1 9 1 6 1 2 bv lPonati et al.l(l2006ah . 
using high-resolution spectropolarimetery. 

A field of order 100 G ha s been reported in the 09.7 su- 
pergiant £ Ori A = HD 37742 dBouret et alJl2008h . but the field 
detection is near the limit of detectability and needs confirma- 
tion. 

IHubrig et alJ (|2008b) carried out the first large survey of O- 
type stars with FORSl, and reported field measur ements of 13 
O star s (including one star previously reported by Hubri g~et alJ 
2007b). Field discoveries were claimed for five of these stars, 
as listed in Table [5] When reduced with our tools, only one 
of these stars still shows a significant magnetic field, the 
06.5f?p star HD 148937. The field of this star has more re- 
cently been detected in three more observations with FORS2 by 
Hubr ig et al.l(l201 lbh . Furthermore, the field has been confirmed 
bv lWade et al.l d201 ll) from observations obtained by the MiMeS 
project. This discovery is clearly real, but it appears that the field 
of this star, always about -250 G, is close to the limit for reliable 
field detection by FORSl in such stars. 

The other O-t ype stars for which fields are reported by 
Hubrig et al. (2008b) appear to be spurious detections. 

Hubr ig et alJ d201 Id) have reported the detection of a mag- 
netic field at just over the 3 cr level in the runaway O-type star 
(, Oph = HD 149757. Our re-reduced measurement of this star 
no longer shows a significant field. This detection is probably 
spurious. 

5.8. X-ray binary star systems 

Under this heading, we consider two different types of possi- 
bly magnetic stars: nearby A-type stars identified with X-ray 
sources in the ROSAT catalogue, and the brilliant and distant 
X-ray source Cyg X- 1 . 

A-type stars are not expected to emit detectable levels of X- 
ray flux. They are not hot enough to have the strong radiatively 
driven winds whose instabilities or collision regions are strong 
X-ray generators, and the A-type stars do not have the deep con- 
vection zones that lead to stellar activity and essentially univer- 
sal coronal X-ray emission. However, about 15% of A-type stars 
are associated with detected X-ray sources. Some of these X-ray 
sources are low-mass, active companions, but it is not clear that 
in all cases this is an adequate explanation. 

ISchroder et al.l (120081) have carried out an extensive survey 
with FORS 1 for magnetic fields in X-ray emitting A-type stars 
in which the companion hypothesis seems doubtful. This sur- 
vey has led to three reported discoveries of fields in A-type 
stars, as listed in Table [5] In the cases of HD 147084 and 
HD 159312, the new reductions show no significant field, and 
these detections are probably spurious. The other detection is in 
HD 148898, which is known to be a magnetic Ap star. Although 
no field has previously b een convincing ly detected in this s tar 
dBorra & Landstreetl 19801) . it is expected dAuriere et a l. 2007) to 
have a weak magnetic field. The field reported by Sch roder et al.1 
(2008) is still present at almost the 3 cr level in our new reduc- 
tions, and may be correctly detected by this observation. 



The system Cyg X-l = HD 226868 is a binary composed of 
a black hole orbiting with an 09.7Iab supergiant, which is los- 
ing mass onto the black hole. FORSl has been used to search 
for the presence of a magnetic fie l d in the opt ical primary star 
of this system bv lKaritskava et al l (120091 1201 Oh . They report 13 
field measurements spread over 1.5 months. Three of their mea- 
surements show apparently significant detections. In our new re- 
ductions, two of these three measurements are still significant at 
about 5 cr level. Two other measurements, not showing a 3 cr de- 
tection in the original reductions, reveal a ^ 3 sigma field detec- 
tion in our new reductions. The field apparently detected in this 
star is very near the reliable detection limit for FORS 1 . We also 
note that due to the target declination, all observations were ob- 
tained with the telescope at a large zenith distance (~ 60°), and 
that in the course of the 1 . 1 h long observations, the instrument 
was rotating by ~ 20°. Should flexures play a role in the mag- 
netic field determinations obtained with FORS, they certainly 
have an impact on these observations. Finally, we must remark 
that we measure an unusually high fraction of circular polariza- 
tion in the continuum (about -0.4 % at 3600 A, linearly increas- 
ing up to about -0.3 % at 5200 A). While this signal is possibly 
real, it could potentia lly be due to cross- talk from linear polariza- 
tion (see Sect. 7.4 of lBagnulo et al.ll2009h . It is clear that further 
observations of Cyg X-l, both in linear and circular polarization, 
would be of considerable interest. 

5.9. Hot subdwarfs 

Detections of kG magnetic fields in two sdO stars, Fei ge 66 
= BP +25° 2534, and BD +75° 325, were reported by Elkin 
(fl996h . 

Magn etic fields were mea sured using FORS 1 in six hot sub- 
dwarfs bv lQ' Toole et al.1 d2005l) . and significant detections were 
claimed for all six stars, as listed in Table [5] Four of these stars 
are B subdwarfs, probably stars that are essentially in the hori- 
zontal branch stage of evolution but have somehow lost virtually 
all their hydrogen envelopes. The other two, CD -31 4800 and 
LSE 153, are O subdwarfs, probably stars on their way from be- 
ing central stars of planetary nebulae to being white dwarfs. 

In our new reductions, the uncertainties of the (B z ) mea- 
surements are slightly reduced compared to the values originally 
published, but no detections at even the 2 cr level are found. This 
survey provides no significant evidence for the presence of mag- 
netic fields in such stars at the level of 1 kG. Furthermore, new 
ESPaDOnS measurements of HP 76431, one of t he sdB stars 
reported to host a kG field bv lQ' Toole et al.l d2005l) . find no in- 
dicati o n of any field ab ove a level of order 100 G dPetit et al.l 
1201 lh . lPetitet all d20 111) have also shown that the FORSl field 
measurement of HP 7643 1 does not support any significant field 
detect ion. Finally, they also report no detection in Feige 66, for 
which Elkin (1996) claimed detection of a kG field. 

IValvavin et al.l J2006) have carried out a search for weak 
(kG) magnetic fields in a number of white dwarfs, and the 
sdO star WP 1036+433, using the spectropolarimeter on the 
Russian 6-m telescope of the Special Astrophysical Observatory. 
Of the three observations of this subdwarf, one reveals a field of 
9.6 ± 2.6 kG, at the 3.7 cr level, which is suggestive of the pres- 
ence of a field, but not yet conclusive. 

At present there is at most marginal evidence for fields of a 
few kG in subdwarfs, and it appears that most of the reported de- 
tections are spurious. However, the total sample studied is very 
small. This is a category of star that would clearly benefit from 
a larger survey with FORS. 
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5.10. Central stars of planetary nebulae 

IJordan et alJ (I2005I) carried out six field measurements of four 
central stars of planetary nebulae (CSPN), and concluded that at 
least two of these stars have highly significant fields in the kG 
range. This result supports the idea that such fields have an im- 
portant shaping effect on the the planetary nebulae themselves. 
It also suggests that important loss of magnetic flux must oc- 
cur between this evolution state and that of white dwarf, since 
white dwarfs with as much magnetic flux as that inferred for the 
CSPN are r elatively rare Thes e results were called into question 
recently bv lLeone et alJ d201 lh . who obtained one new measure- 
ment of each of the CSPNs NGC 1360 and LSS 1362, detecting 
no sig nificant field i n eith er star, and who re-reduced the older 
data of lJordanet al. (2005), also finding no significant fields. 

Our new reductions are complete l y co nsistent with the re- 
reductions carried out by Leone et al. (201 1). None of the orig- 
inal field measurements by Jordan et al. (2005) reveals a signif- 
icant field. As remarked bv lLeone et alJ (1201 lh . there is now no 
significant evidence for coherent magnetic fields in the central 
stars of planetary nebulae, and the current best upper limits on 
(B z ) are at roughly the 1 to 2 kG level. 

5.11. (Weak-field) magnetic white dwarfs 

It has been known for 40 years that strong magnetic fields occur 
in a small fraction of white dwarf stars. The fields can be as large 
as ~ 300 MG, and as small as tens of kG. In recent years, there 
has been considerable interest in establishing how the low field 
end of the distribution of field strengths found in the sample of 
magnetic white dwarfs behaves: do most or all white dwarfs have 
weak (a few kG) fields, or are weak fields as rare as the stronger 
fields are, occurring in only a few percent of white dwarfs? 

A first effor t to a nswer this question was made by 
ISchmidt & Smith! d 1995b . who surveyed nearly 170 DA white 
dwarfs (i.e. white dwarfs with stro ng Balmer line spectra ) using 
the Zeeman effect in the H lines. ISchmidt & Smith! (1 19951) found 
two white dwarfs with fields of order 10 2 kG. More rece ntly, 
a smaller survey was carried out by Fabr ika et alJ (120031) and 
Valv avin et al.l (120061) . using the Russian 6-m telescope, which 
has reported several field detections. The fields detected range in 
size from a few kG to a few tens of kG. 

FOR S1 has made a sub stantial contribution to this search. 
lAznar Cuadrado et al.l J2 004) detected kG fields in three white 
dwarfs. Ijordan et al.l d2007l) "have continued the survey, and iden- 
tified three more stars that they consider probably weakly mag- 
netic (although no single field measurement reached the 3 cr 
level of significance). Our re-re duced data for the three mag- 
netic white dwarfs reported by Aznar Cuadrado et al. (2004) 
still show field values that consistently differ from zero by be- 
tween 3.3 and 6.2 cr. Although none of these measurements is 
completely outside of the range of spurious detections, the fact 
that significant fields are detected in two measurements each of 
WD 0446-789 and WD 2359-434, and that the single measure- 
ment of WD 1 105-048 is significant at the 5.3 cr level, suggests 
that all three field detections may well be correct. However, con- 
firming observations are required. In c ontrast, the ( less signif- 
icant) field measurements reported by IJordan et al l d2007l) all 
have diminished to an insignificant level in the new reductions. 
At present there is no strong evidence for fields in any of these 
stars. 

Note that the sign of the field values has been reversed to 
conform to the convention used in magnetic main sequence stars. 



Since the fields apparently detected with FORS 1 are among 
the very smallest fields known in white dwarfs, it is clear that the 
instrument is a very powerful tool for this kind of work, and that 
further observations would be of considerable value. 

6. Conclusions 

We have investigated the robustness of magnetic field detections 
obtained with the FORS 1 instrument of the ESO VLT. To carry 
out our study, we have developed a semi-automatic data reduc- 
tion procedure and applied it to the entire FORS 1 archive of cir- 
cular spectropolarimetric data. 

We have performed sophisticated Monte Carlo simulations 
(based on the repetition of the entire data reduction process af- 
ter altering the original raw frames to simulate photon noise), 
and concluded that our estimate of the field error bar correctly 
accounts for the photon noise. 

We have performed all our field measurements adopting sev- 
eral different and reasonable data reduction algorithms, and con- 
cluded that the taking certain (somewhat arbitrary) choices dur- 
ing the process of data reduction contributes to the budget of the 
measurement error in a significant way. 

We have determined the magnetic field from the null pro- 
files, and studied its distribution. Our conclusion is that photon- 
noise is not the only source of uncertainties. Among the possible 
causes of non-Poisson noise, we find that small instrument flex- 
ures, negligible in most of the other instrument applications, set 
a lower limit to the accuracy of the field measurements, and may 
be responsible for some of the field detections in the null pro- 
files. 

Overall, non-statistical noise cannot be characterised in de- 
tail, but as a rough guide we estimate that actual error bars may 
be some 50 % higher than estimated from photon-noise, espe- 
cially for high signal-to-noise ratio measurements. 

We have performed a detailed comparison with the results 
previously published in the literature. Even without considering 
the impact of non-photon noise, we are unable to confirm many 
of the previously published new discoveries that were obtained at 
a 3 to 6 cr level, and we have discussed how the incidence of the 
magnetic field in various classes of stars should be revised. The 
fact that our analysis of FORS 1 data does not allow us to confirm 
many of the recent discoveries, leads us to conclude that there is 
no contradiction between the results obtained with FORS 1 and 
those obtained with high-resolution spectropolarimeters. 

FORS is a low resolution spectropolarimeter perfectly suit- 
able for the detection of magnetic fields with ^ 300 G strength. 
Fields weaker than about 100 G are probably out of its reach. 
In general, we argue that any field detection reported at a less 
than 5-6 cr level should be corroborated by repeated observa- 
tions, and, if possible, cross-checked with other instruments. 
While it cannot reach the same accuracy of a high-resolution in- 
strument sitting in a thermally and mechanically isolated bench 
(like ESPaDOnS and HARPSpol), FORS remains an invaluable 
instrument for the detection of magnetic fields in faint and/or 
rapidly rotating stars. 
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Table 5. List of published "marginal" field detections 



Class Star 


MJD 


<B Z > (G) 


Reference 


(B z ) (G) 


Field 


Name 




(published) 




(this work) 


detected ! 


HgMn and PGa stars 












HD 358 


52910.092 


-261±73 


HNS06 


-296+121 


N 




53629.286 


-73+20 


HNS06 


23±33 






53638.205 


-108±23 


HNS 06 


-36±37 




HD 19400 


52852.371 


217+65 


HNS06 


124±85 


N 


HD 65949 


53002.082 


-290±62 


HNS06 


-42+ 79 


N 


HD 65950 


53002.067 


-179+53 


HNS06 


38±74 


N 


HD 175640 


52901.037 


207+ 65 


HNS06 


31 + 64 


N 


Be stars 












HD 56014 


53511.976 


-146±32 


HYP07 


-114+48 


N 


HD 62367 


54549.096 


99±32 


HSSY09 


36+51 


N 


HD 120324 


53869.296 


-80± 24 


HSSY09 


78+111 


N 


HD 148184 


53532.224 


83+21 


HYP07 


113+ 120 


N 




53862.380 


136+16 


HYP07 


-593+522 




HD 181615/6 


53520.415 


38+10 


HYP07 


28+28 


N 




54333.020 


-78±8 


HSSY09 


15+11 






54343 098 


-73+9 


HSSY09 


56+ 10 




HD 209409 


53955.185 


85 ±28 


HSSY09 


-42±46 


N 


HD 224686 


54432.065 


74+24 


HSSY09 




N 


NGC 3766 MG 47 


54549.021 


-234+ 69 


McS08 


-61±60 


N 






-134±42 


T t o o "\ r r\f~\ 

HSSY09 






NGC 3766 MG 200 


54550.375 


-12+42 


McS08 


—1 + 56 


N 






128+40 


HSSY09 






Herbig AeBe stars 












HD 31648 


53296.355 


87+22 


HYS06 


-33±50 


N 


HD 36112 


53331.210 


-149±41 


WBD07 


-73±67 


N 


HD 72 106 A 


53332.296 


195+45 


WDB05 


206+58 


Y 






94+55 


WBD07 






HD 97048 


54609.138 


164+42 


HSSG09 


102+62 


P 


HD 100546 


54610.044 


89+26 


HSSG09 


-43 ±40 


N 


HD 101412 


53062.299 


430+75 


WDB05 


471±136 


Y 






446+ 106 


WBD07 






HD 135344B 


54610.144 


-38+11 


HSSG09 


-14± 17 


N 


HD 139614 


52904.040 


-450+93 


HSY04 


-84±65 


N 






-150+50 


WDB05 










-450+93 


HYS06 










-112+36 


HPY07 








53405.373 


-116+34 


HYS06 


8±25 








-93+ 14 


HPY07 






HD 144432 


53447.352 


-119+38 


HYS06 


-108±22 


P 






-1 1 1± 16 


HPY07 






HD 144668 


52901.007 


166+40 


HPY07 


-90±85 


P 




5461 ?38 


-62+ 18 


HSSfi09 

1 10 0VJU7 


-1 08+ 3? 




HD 150193 


54609.093 


-144+32 


HSSG09 


-238±40 


P 


HD 176386 


54610.272 


— 119+ 33 


HSSG09 


—81 + 41 


N 


HD 190073 


54609.410 


104± 19 


t tci ri s~~* r\r\ 

HSSG09 


-47±57 


Y 


CPD -53 295 


53330.085 


129+ 32 


11 /"Til 1 "X /" \ — T 

WBD07 


132±36 


P 




54610.399 


103+29 


HSSG09 


-5±32 




BF Ori 


53330.277 


-144±21 


WBD07 


-102±45 


N 


p i^ep SlaiS 












HD 16582 


54109.066 


-88+24 


HBD09 


-171±86 


P 




54344.200 


-40± 12 


HBD09 


-1±22 






54344.264 


-44± 12 


HBD09 


47±23 






54345.203 


-42± 12 


HBD09 


-66±24 






54345.246 


-49± 13 


HBD09 


-3±21 






54345.293 


-31+11 


HBD09 


55 ±20 
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Table 5. List of published "marginal" field detections 



Class Star 


MJD 


<B Z >(G) 


Reference 


(B z ) (G) 


Field 


Name 




fmihlishprn 




(\\\\ s worlc^ 

\ LI 11 lj VV V 7 1 IV 1 




HD 50707 


54107.318 


163+52 


HBD09 


54+70 


N 




54345.372 


149+ 19 


HBD09 


67±27 




TTr\ ^ a r~ ■ — 1 r~ 

HD 74575 


54082.341 


i A r\ An 

142+48 


HBD09 


-110+67 


N 




c a 1 r\r\ 1 c r\ 

54109.150 


132+50 


¥ ill r\An 

HBD09 


-316+76 




HD 136504 


c a a a r\r\ci 

54344.998 


-156+34 


t T"n r\AA 

HBD09 


-132+38 


Y 


HD 180642 


C A <~\ A A f\ O A 

54344.084 


166+41 


HBD09 


73+38 


N 


SPB stars 












HD 3379 


53244.402 


272± 57 


HBS06 


84+50 


N 




53245.214 


231 + 47 


HBS06 


-41 + 38 






54344.233 


117+34 


HBD09 


36+36 






54345.189 


155+42 


HBD09 


-14+40 




HD 1 1462 


54344.248 


161 + 46 


HBD09 


23 ±45 


N 


HD 24587 


54086.175 


-353+82 


HBD09 


-140+71 


N 


HD 25558 


54345.264 


105+34 


HBD09 


37+40 


N 


HD 28114 


54106.091 


107±33 


HBD09 


-450+ 104 


N 


HD 28475 


54345.312 


160+48 


HBD09 


7±41 


N 


HD 40494 


54343.426 


94±28 


HBD09 


9±29 


N 


HD 45284 


53252.365 


245+63 


HBS06 


28 ±47 


N 


HD 53921 


52999.137 


-294+63 


HNS06, HBS06 


486±89 


Y 




53630.401 


151 + 29 


HBS06 


153±158 






53631.408 


151 + 21 


HBS06 


154+ 123 




HD 74195 


53138.972 


310+98 


HBS06 


65±83 


P 


HD 74560 


53002.141 


-199+61 


HNS06, HBS06 


259±95 


N 




54108.348 


-198+55 


HBD09 


-68±60 




HD 85953 


53002.152 


-131 + 42 


HNS06, HBS06 


14±55 


N 




54156.096 


97±29 


HBD09 


12±39 




HD 140873 


53151.192 


286+60 


HBS06 


4±56 


N 




54344.011 


99+31 


HBD09 


-29±36 




HD 152511 


54344.116 


649± 43 


HBD09 


595±51 


P 




54608.158 


141 + 26 


HBD09 


66±32 






54609.433 


440± 39 


HBD09 


117±570 






54610.223 


158+28 


HBD09 


118+33 




HD 152635 


54344.041 


-149+36 


HBD09 


-88±30 


N 


HD 160124 


53151.259 


456+60 


HBS06 


-43 ±57 


N 


HD 161783 


53151.281 


376±63 


HBS06 


40±51 


N 




53520.308 


-113+32 


HBS06 






HD 163254 


54344.068 


155+49 


HBD09 


-49±53 


N 


HD 169467 


54345.164 


-182+41 


HBD09 


-74±39 


N 


HD 169820 


53151.312 


239+70 


HBS06 


33±60 


N 


HD 179588 


54343.134 


158+41 


HBD09 


51 + 41 


N 


HD 181558 


53227.184 


236±75 


HBS06 


74±63 


N 




53275.143 


-336+63 


HBS06 


46±55 






54344 1 67 


-104+ 32 


HBD09 


_qn+ 3fi 




HD 183133 


54344.179 


152+38 


HBD09 


81±44 


N 


HD 205879 


54343.226 


150+40 


HBD09 


98 ±49 


N 


HD 208057 


53597.166 


-156+31 


HBS06 


-124±36 


Y 


HD 215573 


52900.080 


165+53 


HNS06, HBS06 


146±70 


N 




Jj Ly 1 .ZZZ 


loU+ J4 


rlDoUu 


A , A A 






O 1 7J.JZ. 1 


_370+90 


HBS06 






Srt stars 
















4/i 1 j 


MluUo 


1 O-i- 1 7 


IN 


Normal B stars 












HD 52089 


54046.339 


-200+48 


HBD09 


-127±60 


p 




54343.389 


-129+34 


HBD09 


-196±37 




6 Car 


54181.025 


266± 52 


HBM08 




N 




54181.035 


-458+118 


HBM08 
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Table 5. List of published "marginal" field detections 



Class Star 


MJD 


<B Z ) (G) 


Reference 


(B z ) (G) 


Field 


Name 




(published) 




(this work) 


detected? 




54181.037 


-394±132 


HBM08 


8+23 






54181.041 


-210+52 


HBM08 








54181.053 


-394+ 125 


HBM08 






NGC 3766 MG 45 


54550.067 


-185+53 


McS08 


-61+54 


P 






-123+40 


HSSY09 






NGC 3766 MG 1 1 1 


54549.020 


54+33 


McS08 


-11 + 40 


N 






112+34 


HSSY09 






NGC 3766 MG 176 


54550.016 


3+31 


McS08 


-4+ 39 


N 






89+28 


HSSY09 






HD 153716 


C A*1 A A f \ C '—I 

54344.057 


124+41 


¥ III r\AA 

HBD09 


-83+50 


N 


HD 179761 


52822.280 


-267+68 


HNS06 


-222+84 


N 


Normal O stars 












HD 36879 


54345.389 


180+52 


HSS08 


33+58 


N 


HD 148937 


54550.416 


-276+88 


HSS08 


-242+ 83 


Y 


iir\ 1 a Ann 

HD 149757 


54609.340 


145+45 


HOS11 


129+58 


N 


HD 152408 


53556.216 


-89+29 


HSS08 


-51+ 114 


N 


lir\ i f f riA/' 

HD 155806 


53531.775 


-115+37 


¥ ¥"* 7 T"W \^ 

HYP07 


70+54 


N 


T IPv 1 A 'I f\ A 

HD 164794 


53594.119 


211 + 57 


HSS08 


169+76 


N 


X-ray stars 












HD 147084 


53975.968 


-48+ 15 


SHS08 


46+25 


N 


HD 148898 


52763.349 


122+29 


SHS08 


184+63 


P 


HD 159312 


53976.178 


241+80 


nTTPAO 

SHS08 


127+78 


N 


HD 226868 


54291.268 


101+ 18 


KBH10 


163+35 


P 




C A *f f A 1 f\ A 

54664.194 


80+23 


KBH10 


260+44 






C A '—1 C\ 1 no 

54678.178 


128+21 


KBH10 


69+ 34 




Hot subdwarfs 












UVO 0512-08 


53058.025 


-1306+161 


OJF05 


155+163 


N 


CPD -64 481 


53058.069 


-885+207 


OJF05 


-72+247 


N 


PG 0909+276 


53058.139 


-1448+222 


OJF05 


-258+220 


N 


CD -31 4800 


53058.215 


-1050+ 161 


OJF05 


-34+ 130 


N 


HD 76431 


53058.255 


-1096+91 


OJF05 


44+ 88 


N 


LSE 153 


53058.347 


-1128+212 


OJF05 


-44+ 156 


N 


Central stars of planetary nebulae 










NGC 1360 


52946.291 


-1343+259 


JWO05 


313+286 


N 




52988.235 


1708+257 


JWO05 


374+262 






52989.060 


2832+269 


JWO05 


564+308 




EGB5 


52988.347 


1992+562 


JWO05 


-1009+773 


N 


LSS 1362 


52989.302 


1891+371 


JWO05 


-H + 394 


N 


Abell 36 


53031.287 


1 \ f l\ , ASS 

1169+466 


JWO05 


1048+379 


N 


Weak-field white dwarfs 












WD 0446-789 


52609.229 


3115+763 


AJN04 


-2765+825 


Y 




52668.087 


/COO 1 , Q-*)Q 

OjZ 1 + yZy 


A T\Tf\/1 

AJINU4 


— JJJJ+ y /4 




WD 1105-048 


52669.305 


-3959+710 


AJN04 


3587+673 


P 


WD 1620-391 


53136.301 


-1580+475 


JAN07 


35+306 


N 


WD 2007-303 


53132.382 


-1093+453 


JAN07 


563+365 


N 


WD 2039-202 


53167.393 


-1297+512 


JAN07 


781 + 598 


N 


WD 2359-434 


52583.025 


-4504+958 


AJN04 


4248+ 878 


Y 




52608.056 


-2919+526 


AJN04 


3211+524 





Key to refer e nce: A JN04: lAznar Cuadrado et all (120041) : JWO 05: iJordan et all (120051) : JAN07 : iJordan et al.1 (l2007|) : HSY04 
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